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geometry,  and  higher  external  impact  indicated  soil-embedded  structure  response, 
which  consisted  of  similar  stages. 

Dimensional  analysis  of  the  model  system  proposes  extrapolation  ratios  and  a 
methodology  to  evaluate  the  distortion  due  to  gravitational  effects,  for  extrapolation 
of  the  small  scale  model  results  into  a  larger  scale  system.  In  cases  of  dynamic 
phenomena  which  are  dominant,  compared  to  the  static  effect  caused  by  gravitational 
acceleration,  and  in  the  lack  of  gravity  scaling,  dimensional  analysis  suggests  the 
distortion  in  the  results  of  the  small  scale  model  are  in  the  range  of  the  experimental 
error,  and  therefore  are  relatively  insignificant. 
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CHAPTER  1 


INTRODUCTION 

LI  Presentation  of  the  Problem  and  Report  Layout. 

This  report  summarizes  a  study  of  the  response  of  a  soil-embedded  structure 
system,  which  is  subjected  to  a  short  impact  applied  at  the  soil  surface.  In  the 
absence  of  a  buried  structure,  those  stresses  induced  in  a  granular  medium  (soil) 
subjected  to  an  external  dynamic  load  (e.g.,  explosive  load)  are  denoted  as  the  "free 
field  stresses".  When  a  structure  is  embedded  in  the  soil  the  field  stresses,  induced 
under  the  same  external  load,  differ  from  the  free  field  stresses,  and  they  now  are 
also  affected  by  the  response  of  the  structure.  The  response  of  the  structure,  which 
is  affected  by  the  stresses  acting  on  it,  is  called  soil-structure  interaction. 

Four  basic  phenomena  were  observed  to  influence  the  total  response  of  the 
system: 

1.  The  external  loading,  which  can  be  either  static  or  dynamic.  In  the  present 
study  the  external  loading  is  an  impact  load  applied  at  the  soil  surface.  That 
generates  a  compression  wave,  which  propagates  into  the  soil  medium.  By 
itself  it  constitutes  the  incident  free  field  stress  distribution. 

2.  A  tension  wave  reflected  from  the  bottom  free  surface  of  the  structural 
element.  In  the  case  of  an  impact  generated  at  the  soil  surface,  that  element 
is  the  roof  of  the  embedded  structure.  Until  the  arrival  of  the  reflected  wave 
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to  the  upper  surface  of  the  roof,  the  incident  wave  is  amplified  at  the  initially 
rigid  soil-structure  interface.  However,  since  the  velocity  of  the  elastic  wave 
in  the  structure  is  relatively  high,  the  reflected  tension  wave  reaches  the  soil- 
structure  interface  at  a  relatively  early  time  of  the  response  and  tends  to 
decrease  the  pressure  level. 

3.  The  structural  element,  which  also  forms  a  physical  boundary  of  the  soil,  gains 
velocity  upon  being  hit.  Thus,  as  long  as  there  is  no  gap  between  the 
structure  and  the  soil  particles,  their  velocity  is  dictated  by  their  boundary, 
that  is,  the  structure.  The  interface  pressure  on  the  structure  is,  therefore, 
influenced  also  by  the  velocity  of  the  interface:  When  the  velocity  of  the  soil 
particles  (at  the  interface)  increases,  the  interface  pressure  decreases,  and  vice 
versa.  The  variation  in  the  interface  pressure  is  considered  relative  to  the 
free  field  pressure  at  the  same  location. 

4.  When  the  external  load  is  static  a  phenomenon  called  arching  has  been 
observed.  Due  mainly  to  its  shear  capacity,  the  soil  transfers  the  load  from 
the  more  flexible  parts  of  the  structure  to  its  more  rigid  parts,  creating  an 
arch-like  stress  field  over  the  structure  that  is  caused  by  the  displacement  of 
the  structure.  Large  scale  experiments  indicate  that  dynamic  soil-structure 
interaction  occurs.  They  also  indicate  the  transfer  of  load  from  the  more 
flexible  parts  of  the  structure  (mid-span  of  the  roof)  to  its  rigid  supports.  This 
behavior  indeed  resembles  arching  as  observed  in  the  static  cases. 
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These  mechanisms  and  the  embedded  structure  behavior  are  coupled  because  they 
influence  the  response  of  the  structure,  and  are  also  initiated  by  it.  Their  influence 
may  differ  according  to  their  timing.  For  example,  prior  to  the  development  of  a 
structural  displacement  sufficient  to  cause  arching,  the  structure  will  for  dynamic 
loads,  accelerate  and  gain  velocity,  thereby  defining  a  moving  boundary  for  the 
pressure  wave  propagating  in  the  soil.  The  possibility  of  initiating  the  development 
of  arching  may  occur  as  a  result  of  this  wave.  At  some  subsequent  time  the  structure 
might  rebound,  gaining  an  upward  velocity,  and  altering  the  stress  field  which  was 
initially  created  by  the  prior  occurrence  of  the  mechanisms  described  above. 

In  the  present  study,  better  understanding  of  these  fundamental  phenomena  is 
sought  through  a  small  scale  experimental  system.  By  closely  monitoring  the  loads 
over  the  buried  structure  and  recording  its  associated  velocity  and  displacement, 
appropriate  conclusions  may  be  drawn  regarding  the  existence  and  influence  of  these 
mechanisms  on  the  soil-structure  system. 

The  second  part  of  this  chapter  is  a  literature  survey  detailing  the  various 

approaches  to  the  problem  and  the  pertinent  observations  and  conclusions  drawn  by 

other  researchers.  The  experimental  part  of  this  research  is  described  in  Chapter  2, 
* 

which  includes  a  description  of  the  experimental  setup,  the  test  plan,  results  and  the 
conclusions.  The  last  section  of  Chapter  2  discusses  the  value  of  the  results  obtained 
in  the  small  scale  system,  as  reflected  by  the  pertinent  scaling  considerations. 
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Based  on  the  experimental  results,  a  non-linear,  dynamic,  finite  element  computer 
code,  "SAMSON2",  was  calibrated  and  used  for  further  numerical  analysis.  The  code, 
the  numerical  parameters  used,  and  the  calibration  process  are  described  in  Chapter 
3.  The  parameters  of  the  materials  models  were  calibrated  according  to  the 
experimental  results  of  the  small  scale  system.  Further  output  of  the  numerical 
analysis  of  the  system,  such  as  the  shear  stresses  in  the  soil,  or  the  vertical  stresses 
in  the  soil  adjacent  to  the  structure,  was  obtained.  These  records  were  not  measured 
in  the  experiments,  and  contributed  additional,  and  better  insight  into  the  study  of 
the  soil-structure  system. 

Chapter  4  presents  the  application  of  the  computer  program  to  further  analytical 
study  of  the  similar,  but  larger  scale  problem  of  a  soil-embedded  structure  system 
subjected  to  an  impact  load.  A  reference  field  test  which  had  planar  geometry  of 
larger  dimensions,  and  a  higher  external  load,  generated  by  an  explosive  charge,  was 
simulated  numerically.  The  significant  difference  between  this  soil-structure  system 
and  the  smaller  one  was  the  level  of  the  external  load.  This  difference  was  reflected 
in  the  results  of  the  analysis  by  the  higher  initial  loading  of  the  structure.  However, 
the  mechanisms  which  occurred  in  the  smaller  scale  system,  were  also  observed  in 
the  larger  reference  system. 

The  small  scale  laboratory  experiments  together  with  the  subsequent  numerical 
analysis  indicate  patterns  of  behavior,  which  are  typical  to  the  response  of  a  shallow 
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embedded  structure,  subjected  to  an  external  impact.  They  emphasize  the  dominant 
effect  of  the  early-time  response,  and  the  succeeding  influence  of  the  later 
mechanisms.  These  are  summarized  in  the  closure.  Chapter  5. 
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L 2  Literature  Survey 

The  response  of  a  structure  to  external  load  is  measured  by  its  deflections,  strains, 
and  ultimately,  by  its  capacity  to  sustain  the  load.  Embedded  structures  are  usually 
constructed  inside  an  excavation  in  the  soil,  which  is  then  backfilled  around  and  over 
the  structure.  In  this  case  the  load  is  usually  not  acting  directly  on  the  structure,  but 
on  and  through  the  soil  in  which  it  is  embedded.  Load  may  be  applied  to  the  soil 
surface,  as  in  the  case  of  a  railroad  over  a  conduit,  or  it  may  be  an  external  explosive 
load.  It  can  also  be  a  lateral  load  which  is  typical  for  earthquakes.  The  response  of 
the  embedded  structure  is  caused  and  affected  by  the  load  induced  on  it.  However, 
because  that  load  is  applied  through  the  soil  it  is  also  affected  by  the  response  of  the 
soil  and  by  its  strain  field,  and  the  state  of  the  soil  is  influenced  by  the  response  of 
the  structure  embedded  in  it.  This  coupled  behavior  is  called  soil-structure 
interaction. 

Early  works  on  soil-structure  interaction  related  to  embedded  structures  dealt  only 
with  external  static  loads.  In  his  book,  Terzaghi(471  analyzed  experimental 
observations  which  indicated  the  formation  of  failure  shear  planes  above  a  horizontal 
trap-door  that  was  displaced  uniformly  downward  (Fig.  1-1).  This  example 
demonstrated  the  influence  of  the  structural  displacement  on  the  load  acting  on  it: 
Assuming  development  of  shear  stresses  along  the  failure  shear  planes,  and  their 
relation  with  the  vertical  and  horizontal  stresses,  Terzaghi  suggested  the  following 


expression  for  the  average  vertical  pressure,  p,  over  the  deflecting  trap-door: 
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where:  p  is  the  average  pressure  over  the  trap-door  (p  =  av(Z=D)). 

B  is  half  of  the  trap-door  span. 

Y  is  the  soil  weight  density. 

C  is  the  soil  cohesion. 

K  is  the  lateral  pressure  coefficient. 

^  is  the  internal  friction  angle. 

D  is  the  depth  from  the  surface  to  the  structure, 
q  is  a  uniform  external  pressure  on  the  surface. 

Although  equation  1.1  was  developed  for  a  particular  case  and  under  certain 
assumptions,  it  is  representative  of  the  influence  of  a  buried  structure’s  deflection  on 
the  average  load  acting  on  it: 

a.  The  larger  the  depth  to  span  ratio  (D/B)  the  smaller  the  load. 

b.  Beyond  a  certain  depth  to  span  ratio  (-5)  the  influence  of  the  external  load, 

q,  disappears,  and  the  average  load  becomes  constant  and  equal  to  . 

JTtan<t> 
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c.  The  average  load  over  the  structure  decreases  as  the  cohesion  and  angle  of 
internal  friction  of  the  soil  increase. 

These  conclusions  remain  essentially  the  same  when  the  expression  for  the  average 
pressure,  p,  is  compared  with  the  free  field  pressure  at  the  same  level,  q+ yD.  Such 
a  comparison  shows  that  the  pressure  level  over  the  relatively  flexible  structure  to  be 
less  than  that  of  the  free  field.  To  maintain  vertical  equilibrium,  the  pressure  at  both 
sides  of  the  structure  must  increase. 

Using  a  similar  approach11,2’34’411  exponential  decrease  of  the  pressure  over  a  buried 
structure  has  also  been  suggested  by  other  researchers.  This  approach,  which 
assumes  that  shear  planes  have  already  developed  in  the  soil  provides  a  basic 
understanding  for  the  mechanism  causing  the  load  reduction  above  a  buried  structure 
which  settles  relatively  to  the  surrounding  soil,  and  give  an  estimate  of  the  overall 
load  acting  on  the  structure.  The  failure  mechanism  is  one  in  which  shear  stresses 
act  along  the  shear  planes  thereby  reducing  the  load  acting  on  the  structure  below. 

Formation  of  failure  shear  planes,  however,  is  not  essential  for  the  existence  of 
this  phenomenon  of  pressure  decrease  over  the  structure.  As  it  was  observed  in 
experimental  studies  with  static  load121’23’33’26’27,451,  the  pressure  acting  on  a  relatively 
flexible  structure  decreases  while  the  pressure  at  both  sides  increases  even  before  the 
soil  reaches  a  state  of  failure.  Numerical  analyses114’20,44,551  show  that  the  distribution 
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of  the  principal  compression  stresses  in  the  soil  above  the  structure  constructs  an 
arch  action  (Figure  1-2).  Thus,  this  phenomenon  of  decreasing  pressure  over  a 
relatively  flexible  buried  structure  is  commonly  referred  to  as  positive  or  active 
arching.  In  the  opposite  case  when  the  structure  is  relatively  rigid,  the  pressure  over 
the  structure  (or  rigid  structural  elements)  increases  (relative  to  the  adjacent 
medium).  This  case  is  referred  to  as  negative  or  passive  arching.  It  was  further 
found,  in  the  studies  mentioned  above,  that  arching  in  soils  is  enabled  mainly  through 
shear  stresses,  which  at  a  high  level  of  load  and  structural  deflection,  generate  failure 
shear  planes  in  the  soil.  Increasing  depth  of  burial  (DOB)  causes  increasing  arching, 
i.e.,  decrease  of  the  load  over  the  flexible  structure,  yet  arching  was  also  observed  at 
a  DOB  as  shallow  as  20%  of  the  structure  span,  L1191. 

The  extensive  research  of  the  static  phenomenon  lead  to  its  fundamental 
understanding  and  has  been  used  as  a  basis  for  the  research  of  the  behavior  of  a  soil- 
embedded  structure  subjected  to  a  dynamic  external  load.  Kiger  and  Getchell 
performed  a  series  of  tests  in  which  reinforced  concrete  box-type  structures  where 
buried  under  a  relatively  shallow  DOB  (<L/2)  of  two  different  backfills,  and 
subjected  to  external  explosive  loads  generated  by  High  Explosive  Simulation 
Technique  (HEST)[18,19’25J.  These  notable  experiments  contributed  important 
information  to  the  research  of  the  behavior  of  shallow  buried  structures.  In  three 
tests  (HEST  1,  3,  and  4)  the  roof  of  the  structure  was  damaged  but  did  not  totally 
collapse.  In  these  tests  the  pressure  records  that  were  measured  at  the  roof  center 
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were  different  from  those  that  were  measured  over  the  supporting  walls.  The 
pressure  over  the  center  decreased  and  was  lower  than  that  over  the  walls.  This 
difference  occurred  late  in  the  response,  within  1  to  5  msec  after  the  roof  was  hit  by 
the  compression  wave  and  was  reflected  in  a  similar  way  on  the  corresponding 
calculated  impulse  records.  These  results  were  obtained  at  a  DOB  of  L/2,  but  even 
when  the  DOB  was  as  shallow  as  02L  (HEST  4)  the  impulse  over  the  roof  center 
was  reported  to  be  less  than  one  third  of  the  impulse  over  the  walls.  In  the  former 
case  the  response  of  the  roof  was  flexural,  while  a  mixed  mode  (flexural  and  shear) 
response  was  reported  in  the  latter. 

In  the  only  two  tests  where  the  roof  totally  collapsed  the  external  peak  pressure 
was  well  above  that  of  the  other  tests  and  failure  occurred  at  a  time  much  shorter 
than  the  longest  calculated  period  of  the  roof.  The  failure  mode  was  dominated  by 
strong  shear  near  the  walls.  It  was  believed  that  in  these  tests  early  failure  occurred 
under  the  high  level  of  loading,  before  any  interaction  between  the  structure  and  the 
soil  could  become  significant  enough  to  prevent  it  or  to  reduce  its  severity.  In  the 
other  cases  the  authors  considered  relief  waves  reflecting  back  from  the  concrete-air 
interface  on  the  bottom  face  of  the  roof,  as  a  possible  cause  for  the  fast  reduction 
of  the  interface  pressure  at  the  roof  center.  They  dismissed  it,  however,  and 
attributed  the  explanation  to  dynamic  soil  arching  as,  indeed,  transfer  of  pressure 
from  the  flexible  part  of  the  structure  to  its  more  rigid  elements  is  indicative  of 
arching. 
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A  late  time  (relative  to  the  initial  structural  response)  effect,  which  was  attributed 
to  dynamic  soil  arching,  was  reported  by  Kiger,  Eagles,  and  Baylot1271.  The  dynamic 
part  of  their  experimental  program  consisted  of  a  reinforced  concrete  slab  embedded 
in  the  soil  with  a  depth  to  span  ratio  of  0.5,  and  subjected  to  a  surface  Foam  HEST 
overpressure.  They  measured  a  relatively  higher  impulse  over  the  embedded  slab 
supports  that  was  interpreted  as  possible  soil  arching.  Similar  experimental 
observation  is  pointed  out  in  tests  that  were  made  with  buried  reinforced  concrete 
arches  subjected  to  explosive  overpressure1301.  Late  time  behavior  of  buried  arch-type 
structures  v/as  found  to  be  affected  by  soil-structure  interaction,  especially  in  inducing 
relative  motions  between  the  arch  walls  and  the  floor  after  damage  initiates  at  the 
arch-to- floor  connection.  It  was  also  found  that  the  properties  of  the  soil-structure 
interface  can  have  an  influence  on  the  structural  performance  by  inducing  friction 
loads. 

Windham^541  used  the  finite  element  method  to  investigate  the  dynamic  response 
of  a  box-type,  40  ft  by  25  ft  (122  m  by  7.6  m)  structure.  He  examined  the  influence 
of  the  backfill  parameters  on  the  behavior  of  the  system.  The  late-time  pressure 
distribution  over  the  roof  showed  a  relatively  uniform  distribution  over  the  center 
portion  of  the  roof,  and  a  sharp  normal  stress  concentration  over  the  walls.  No 
significant  difference  between  dense  and  loose  backfill,  of  the  level  or  distribution 
of  the  pressure,  is  reported  in  this  work.  The  only  noticeable  influence  of  the  backfill 
density  is  on  the  normal  stresses  on  the  regions  of  the  roof  located  over  the  walls  and 


12 


on  the  side  walls  at  the  roof  level,  which  were  higher  under  the  loose  backfill. 

Soil  backfill  above  the  structure  affects  its  frequency  response  through  stiffening, 
and  by  an  added  mass  effect.  At  a  shallow  DOB  of  about  20%  of  the  roof  span  the 
added  mass  of  the  soil  appears  to  be  dominant  in  its  effect  on  the  response  of  the 
structure,  as  shown  and  discussed  by  Dallriva  and  Kiger1131,  who  also  observed  in  their 
experiments  an  increase  in  the  system’s  damping  when  soil  is  added  as  backfill. 

It  should  be  noted  that  the  decrease  of  pressure  which  was  observed  in  the 
dynamic  experiments  at  a  very  shallow  DOB,  and  with  a  clay  backfill,  can  also  occur 
due  to  the  velocity  difference  between  the  structure  and  the  medium,  at  least  at  the 
beginning  of  the  roof  motion.  This  interface  pressure  is  governed  by  the  relative 
velocity  between  the  structure  and  the  medium  in  which  it  is  embedded,  and  by  the 
acoustic  impedance ,  pc  of  the  medium.  The  acoustic  impedance  gives  a  relation 
between  the  pressure  induced  in  a  medium  and  its  particle  velocity,  v.  For  example, 
the  relation  p=(pc)v  can  be  applied,  for  a  plane  wave,  at  the  boundary  of  the 
medium,  such  as  the  soil-structure  interface. 

This  approach  was  used  by  Weidlinger  and  Hinmant51J,  who  employed  this 
pressure-velocity  relation  with  a  nonlinear  response  of  the  structure  into  a  single 
degree  of  freedom  system  (SDOF).  They  also  allowed  for  the  formation  of  a  gap 
between  the  sand  particles  and  the  structure  by  setting  the  interface  pressure  to  zero 
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when  such  a  gap  occurs,  and  suggested  that  the  pressure  relief  over  the  roof  center 
might  be  denoted  as  dynamic  arching.  This  notation,  however,  should  be 
distinguished  from  static  arching  which  is  driven  by  a  different  mechanism. 

With  the  exception  of  finite  element  computations,  numerical  methods  which  deal 
with  the  dynamic  response  of  buried  structures,  usually  consider  the  inertial  effect 
and  the  conditions  at  the  soil-structure  interface,  or  the  influence  of  the  acoustic 
impedance,  but  do  not  include  in  their  consideration  the  soil  arching  effect116,221. 
Nevertheless,  these  papers  usually  give  good  agreement  of  their  results  with 
experimental  or  other  numerical  results,  thus  questioning  the  feasibility  of  soil 
arching,  or  the  level  of  its  effect,  on  the  dynamic  behavior  of  a  soil-embedded 


structure  systems. 


Fig.  1-2.  Pattern  of  the  principal  stresses  for  DOB/span=l  in  Yankelevsky’s  numerical  analysis^. 


CHAPTER  2 


EXPERIMENTS 

2.1  Description  of  the  Experimental  System. 

Full  or  large  scale  experimental  research  of  the  response  of  buried  structures  to 
external  impact,  has  the  advantage  of  employing  realistic  parameters:  The  external 
impact  can  be  generated  by  actual  blast  load,  local  soil  may  be  used  for  the  backfill, 
and  the  size  of  the  structure  is  either  full  scale  or  close  to  it.  The  difficulties  that  are 
involved  in  such  programs  are  their  relative  high  cost,  and  inability  to  precisely 
determine  those  parameters.  For  example,  it  is  difficult  to  accurately  measure  the 
pressure  which  is  generated  by  a  blast  or  to  delineate  its  distribution. 

The  purpose  of  the  current  research  was  to  focus  on  the  response  of  the  structure, 
and  to  study  the  basic  mechanisms  that  are  involved  in  a  soil-embedded  structure 
system  subjected  to  external  impact  by  eliminating  the  number  of  its  unknown 
parameters.  The  approach  that  was  selected  was  to  employ  an  experimental  program 
which  consisted  of  small  scale  soil-structure  system  (Figure  2-1).  Such  a  system 
enabled,  within  a  relatively  low  budget,  detailed  measurements  and  information 
which  provide  an  insight  to  the  problem,  and  are  complementary  to  similar  larger 
scale  experiments. 
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2.1.1  Experimental  Setup. 

The  experimental  setup  was  developed  and  constructed  in  the  research  laboratory 
of  Northwestern  University181.  A  circular  steel  tank,  42  in  (107  cm)  in  diameter  and 
about  the  same  height,  contained  20-30  dry  Ottawa  sand.  A  steel  pipe  frame  was 
mounted  around  and  above  the  tank  and  supported  an  electromagnet  which  was 
placed  along  the  height  of  the  frame  at  variable  altitudes.  A  4  pound  (17.8  N),  3  in 
(7.62  cm)  diameter,  steel  ball  was  used  as  an  impactor.  By  dropping  the  ball  on  a 
relatively  rigid  target  plate,  resting  on  the  soil  surface,  a  short  impact  was  generated, 
which  then  propagated  into  the  soil.  The  minimal  height  of  drop  was  25  in  (63.5  cm) 
and  the  maximum  was  90  in  (228.6  cm).  The  height  of  drop  was  measured  from  the 
target  plate  to  the  bottom  of  the  impactor  ball.  The  target  plate  was  made  of 
aluminum,  0.5  in  (1.24  cm)  thick  and  12  in  (30.48  cm)  in  diameter,  sufficiently  rigid 
to  distribute  the  concentrated  impact  generated  by  the  steel  ball  into  the  soil,  and 
light  enough  to  avoid  excessive  artificial  static  load.  Extensive  investigation  of  the 
target  plate  response  was  done  by  Chen1®1. 

The  experimental  system  had  an  axially  symmetric  configuration  that  enabled  a 
convenient  experimental  as  well  as  a  later  numerical  analysis.  Accordingly,  the 
embedded  small  scale  structure  was  a  vertical  cylinder,  which  had  an  external 
diameter  and  a  height  of  6  in  (15.24  cm).  This  experimental  research  was  focused 
on  the  response  of  the  roof  of  the  structure.  The  circular  roof  slab  was  supported 
by  a  plexiglas  cylindrical  wall  with  a  circular,  1  in  (2.54  cm)  thick,  floor.  In  the  first 
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stage  of  the  experimentation  a  linear  elastic  plexiglas  slab  was  used  for  the  roof.  It 
was  0.23  in  (0.58  cm)  thick,  and  was  simply  supported  by  a  0.25  in  (0.64  cm)  thick 
wall.  In  the  second  and  main  phase  of  the  experimentation  the  roof  slab  was  made 
of  micro-concrete  and  reinforced  with  micro-reinforcement.  The  preparation  and 
properties  of  these  slabs  are  described  later  in  the  text,  together  with  the  results. 
The  micro-concrete  roof  slabs  were  supported  by  a  0.5  in  (1.27  cm)  thick  wall,  which 
together  with  the  1  in  (2.54  cm)  thick  floor  formed  a  relatively  rigid  structure,  and 
enabled  the  structural  behavior  to  be  focused  on,  and  dominated  by  the  response  of 
the  roof. 

2 .12  Equipment. 

Measurements  of  load  and  motion  were  performed.  It  is  believed  that  by 
obtaining  and  analyzing  the  pressure  distribution  and  the  related  structural  motion 
(velocity  and  displacement),  the  mechanisms  involved  in  this  soil-structure  system  can 
be  studied  and  understood.  Modified  ELF-500-100  "Entran"  load  cells  were  used  to 
measure  the  pressure  at  the  points  where  they  were  located.  Interpolation  of  the 
pressure  records  acquired  by  load  cells  at  different  locations  provided  the  pressure 
distribution.  Description  of  the  load  cells  and  the  modification  needed  in  order  to 
use  them  in  this  particular  system,  is  given  in  Appendix  A. 

The  motion  of  the  roof  center  was  recorded  by  a  Linear  Variable  Differential 
Transformer  (LVDT).  This  transducer  demonstrated  reliable  displacement  output 
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as  shown  and  discussed  in  Appendix  B.  The  time  derivative  of  the  LVDT  output 
signal  provided  the  velocity  record.  It  should  be  noted  however,  that  the  LVDT 
measured  the  motion  of  the  roof  relative  to  the  structure  floor.  In  addition  to  the 
LVDT  two  4374  "Bruel  &  Kjaer"  accelerometers  were  mounted  at  the  roof  center 
and  on  the  floor.  The  lower  accelerometer  provided  the  floor  motion  record,  which 
was  also  assumed  to  be  the  motion  of  the  whole  structure,  while  the  upper 
accelerometer,  after  double  integration,  provided  the  absolute  displacement  record 
at  that  point.  The  signal  from  the  accelerometer  was  used  to  indicate  the  exact 
starting  time  of  the  roof  motion,  as  its  initial  response  was  faster  than  that  of  the 
LVDT.  Because  the  double  integrated  acceleration  signals  were  less  accurate  than 
the  LVDT  signal,  they  were  considered  to  be  more  indicative  rather  than  precise 
records.  This  feature  is  also  discussed  in  the  plexiglas  slab  test  results  (see  2.3.2)  and 
demonstrated  in  Appendix  B. 

Acquisition  of  the  output  signals  from  the  transducers  was  performed  with  a  4094 
Nicolet  Digital  Oscilloscope,  and  two  2232  Tektronix  Digital  Storage  Oscilloscopes. 


20 


2*2  Test  Elan* 

The  soil-structure  experimental  program  included  four  main  stages: 

a.  Experiments  intended  to  measure  the  external  impact  input  generated  under 
the  target  plate  by  the  drop  of  the  steel  impactor  ball. 

b.  Free  field  measurements. 

c.  Experiment  with  a  plexiglas  roof  slab. 

d.  Experiments  with  reinforced  micro-concrete  roof  slabs. 

22.1  External  Impact  Measurements. 

The  external  load  that  developed  under  the  target  plate,  as  well  as  the  free  field 
stress  distribution,  are  two  important  reference  records  for  this  research.  The 
external  load,  i.e.  the  pressure  generated  beneath  the  target  plate,  is  the  input  to  the 
soil-structure  system,  and  once  determined  experimentally,  was  used  as  an  input  for 
the  subsequent  numerical  analysis.  It  was  also  used  as  a  reference  for  the 
experiments  that  were  performed  during  succeeding  stages,  which  involved  the 
embedded  structure.  In  these  experiments  a  measurement  of  the  pressure  generated 
under  the  target  plate  center  was  taken  to  confirm  that  the  external  input  was  similar 
to  that  of  the  free  field  under  similar  drop  of  height  of  the  impactor.  This  record, 
together  with  the  free  field  records,  well  define  the  experimental  system. 


With  no  buried  structure,  the  target  plate  and  four  load  cells  distributed  beneath 
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it  were  placed  on  the  sand  surface,  Figure  2-2.  The  normal  pressure  was  recorded 
at  the  center  (load  cell  T),  1.75  in  (4.45  cm)  away  from  the  center  (load  cell  ’A’), 
3.25  in  (826  cm)  away  from  the  center  (load  cell  ’B’),  and  5  in  (12.7  cm)  away  from 
the  center  (load  cell  ’C’).  With  the  assumption  of  axisymmetry  a  fourth  order 
polynomial  was  used  for  the  radial  load  distribution  under  the  target  plate,  as 
follows: 

p(r,t)=  a4r4  +  +  ajr2  +  axr+  a^  (2.1) 

where  the  requirement  for  axisymmetry,  yields  : 


djl 

dr 


=  at  =  0 

r*0 


(2.2) 


Hence:  p(r,t)  =  a4r4  +  a3r3  +  fy*2  +  &o 


(2.3) 


By  substituting  the  recorded  data  p(0cm,t),  p(4.45cm,t),  p(8.26cm,t),  and  p(12.7cm,t), 
the  following  matrix  is  obtained: 


0  0  0  1 

a4 

9.38  536  3.06  1 

^3 

111.57  34.33  10.56  1 

a 2 

625.00  125.00  25.00  1 

(2.4) 
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from  which  a4,  a3,  aj  and  ao  are  determined  as  follows  : 


a4 

-0.038  0.067  -0.036  0.007 

4 

= 

0.344  -0.553  0.243  -0.035 

a. 

-0.813  1.088  -0.316  0.040 

I 

3o 

1.000  0.000  0.000  0.000 

pfl.^",t 
p(3.25" 
P(5",t) 


::D 


(2.5) 


It  should  be  noted  that  a4,  a3,  %  and  a*,  are  also  functions  of  time.  Thus,  based  on 
the  recorded  vertical  pressure  (and  assuming  a  fourth  order  distribution,  according 


to  equation  2.3),  the  full  spatial  and  time  distributions  of  the  external  load  are 


obtained  as: 


a4  =  a4ft) 
a3  =  a3(t) 
a*  =  a^t) 
3o  —  ao(0 


P(M) 


Figures  2-3-a  and  2-3-b  show  this  distribution  under  a  height  of  drop,  H  =  25  in 
(63.5  cm).  This  pressure  distribution  was  used  later  as  an  input  for  the  finite  element 
numerical  analysis,  thereby  eliminating  the  need  for  numerical  modeling  of  the  load 
input  transmitted  from  the  impactor  ball  at  the  center  of  the  target  plate  and 
decreasing  the  error  that  might  have  been  introduced  by  doing  so.  The  duration  of 
the  total  loading  pressure  wave  as  seen  in  figure  2-3-a,  generated  under  the  target 
plate,  was  about  1  msec.  The  pressure  under  the  plate  center,  however,  had  a 
duration  of  less  than  0.5  msec,  and  a  peak  of  about  325  psi  (2.2  MPa).  It  was  the 
vibration  of  the  plate  that  transmitted  the  pressure  wave  from  the  center  to  its  edge, 
where  it  had  a  peak  of  about  100  psi  (0.7  MPa)  at  about  0.6  msec  (Fig.  2-3-b). 
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2 .2.1.1  Approximation  for  the  Total  Load  Induced  under  the  Target  Plate. 

The  purpose  of  the  following  discussion  is  to  verify  the  external  input  load 
measurements  by  approximating  the  total  load  induced  under  the  target  plate.  When 
dropped  from  a  height  H,  the  impactor  ball  hits  the  target  plate  with  a  velocity  equal 
to  (2gH)1/2,  where  g  is  the  acceleration  of  gravity.  The  softer  aluminum  plate 
undergoes  some  plastic  deformation.  The  contact  area,  as  measured  on  the  target 
plate,  had  an  average  diameter  of  3/16  in  (0.48  cm).  The  yield  stress  of  aluminum 
is  8000  psi  (55.2  MPa),  therefore  a  force  of  about  200  pounds  (890  N)  was  applied 
to  cause  local  yielding  and  plastic  deformation  at  the  contact  area.  The  ball  and  the 
plate  move  together  for  a  period  of  time  tc,  which  was  measured  by  Chen,8).  For 
H=25  in  (63.5  cm)  the  contact  time,  t^  was  equal  to  0.77  msec.  Then,  the  ball  was 
observed  to  rebound  off  the  plate.  Assuming  that  the  ball  and  the  plate  move 
together  after  the  local  yielding  occurs,  by  momentum  conservation  their  initial 
velocity,  V0  is: 

V0  -  (2.6) 


where  Mb  and  Mp  are  the  masses  of  the  ball  and  of  the  plate,  respectively.  The 
force  induced  under  the  plate  is  assumed  to  be  generated  by  the  weight  of  the  ball 
and  of  the  target  plate,  and  by  the  resistance  of  the  sand  to  their  deceleration. 
Consider  an  equivalent  SDOF  system  with  the  following  parameters  (figure  2-4): 
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Lumped  mass,  M  =  Mb  +  Mp 
Equivalent  resistance  of  the  sand,  K. 

Although  it  is  evident  that  the  response  of  such  system  is  not  linear  and  that  it  is 
damped,  because  of  the  short  duration  of  the  ball-plate  contact  (less  than  1  msec) 
a  linear  undamped  analysis  may  be  sufficient  for  such  an  estimation.  The  equation 
of  motion  is: 

MX  +  Kx  =  Mg  (2-7) 


where  x  =  x(t)  is  the  vertical  displacement.  The  initial  conditions  are: 


m 

m 


K 


Vo 


Mb+Mp 


yfigH 


(2.8) 


The  solution  for  the  displacement  x(t)  is: 

x(t )  =  — sintor  +  — — (1  -cosujr)  +  (2.9) 

o )  K  K 


where  w2=— = — — — 
M  Mb+Mp 
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According  to  this  model,  the  total  load  applied  under  the  target  plate,  FMnd  is: 

KV 

F.  =  Kx  =  W+  WJl-cosot)  +  — -sinor  (21°) 


where  W  denotes  weight  (of  the  ball  or  of  the  plate). 

The  equivalent  parameter  K  can  be  determined  semi-empirically,  since  the 
duration  of  the  load  under  the  target  plate,  was  measured  and  recorded,  and 
Ftand(t  =  td)  can  therefore  be  set  to  equal  Wp,  the  static  load  of  the  plate.  Note  that 
according  to  equation  2.10,  Wp  is  also  the  initial  load,  Fsand(t  =  0)=Wp.  This 
additional  condition  allows  the  semi-empirical  determination  of  K.  Once  K  is 
obtained  it  may  be  substituted  into  equation  2.10  to  obtain  F*and(t)  for  t<td. 

The  loading  duration,  td  for  a  height  of  drop  H=25  in  was  0.4  to  0.5  msec  (figure 
2-3-a  and  2-3-b).  Note  that  it  was  of  the  same  order  of  magnitude  as  the  contact 
time,  measured  by  Chen1®1.  The  weight  of  the  target  plate,  Wp,  and  of  the  ball,  Wb, 
was  5.4  pounds  and  4  pounds,  respectively.  For  td  =  0.4  msec  equation  2.10  yields 
K=  1501000  lb/in  and  for  td=0.5  msec  K= 961000  Ib/in.  The  measured  and 
interpolated  pressure  distributed  under  the  target  plate  was  integrated  numerically 
over  its  area  to  produce  the  total  load.  The  experimental  value  is  plotted  together 
with  the  analytical  prediction  (for  td=0.4  and  0.5  msec)  in  figure  2-5.  As  it  can  be 
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seen  the  SDOF  model  gives  a  reasonable  theoretical  prediction  for  the  total  load 
exerted  under  the  target  plate.  It  should  also  be  noted  that  the  experimental  value 
was  obtained  from  load  cells  which  were  located  at  a  maximum  distance  of  5  in  from 
the  plate  center.  Since  no  extrapolation  was  done  (see  previous  section),  the 
experimental  total  load  does  not  reflect  the  load  contributed  from  the  external 
perimeter  of  the  plate  (which  had  a  radius  of  6  in). 

This  result  is  another  verification  for  the  pressure  measurements  as  recorded  by 
the  setup  of  load-cells  described  in  Appendix  A. 

2 22  Free  Field  Measurements. 

The  vertical  free  field  stress  distribution  under  the  target  plate  is  important  to 
provide  a  reference  for  comparison  to  the  distribution  of  the  stresses  when  a  buried 
structure  is  present.  Furthermore,  the  free  field  distribution  provides  a  reference  to 
check  the  numerical  model  used  for  the  soil  constitutive  relationship  in  the  finite 
element  analysis. 

With  no  structure  present,  load  cells  buried  in  the  soil  measured  the  vertical 
pressure  that  developed  under  external  impact  similar  to  the  impact  generated  later 
in  the  following  stages  of  the  experimental  program.  Figure  2-6  describes  the  test 
configurations  which  were  used  to  measure  the  free  field  vertical  stresses  2.5  in  (6.35 
cm)  and  3  in  (7.62  cm)  beneath  the  target  plate  at  three  locations:  directly  under  the 
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center,  2  in  (5.08  cm)  and  2.5  in  (6.35  cm)  away  from  the  center.  Figures  2-7-a  and 
2-7-b  show  some  results  obtained  from  these  tests,  which  together  with  the  external 
load  measured  under  the  target  plate,  complete  the  reference  data.  Additional 
pressure  records  that  were  obtained  in  these  tests  are  shown  together  with  the  soil- 
structure  results  from  the  next  two  stages  of  the  experimental  program. 

2 23  Experiment  with  an  Elastic  Plexiglas  Roof  Slab. 

Once  the  input  loads  and  the  free  field  stress  distribution  have  been  established, 
the  soil-structure  system  consisting  of  an  elastic  buried  structure  was  tested  (Figure 
2-8).  This  experiment  was  performed  with  a  plexiglas  roof  in  order  to  obtain  a  full 
response  of  the  structure  in  its  linear  elastic  range.  It  was  also  used  for  the 
inspection  of  the  setup  and  its  measuring  and  recording  equipment. 

The  plexiglas  cylinder  caped  with  a  plexiglas  roof  slab,  and  fully  instrumented  with 
the  transducers  described  above,  was  placed  in  the  soil  under  the  target  plate.  The 
plexiglas  slab  was  simply  supported,  and  since  there  was  no  attempt  to  reach  its 
failure,  the  external  impact  was  generated  by  the  lowest  height  of  drop  of  the 
impactor,  25  in  (63.5  cm).  The  slab  thickness  was  0.23  in  (0.58  cm),  and  the  depth 
of  sand  cover  above  it,  and  below  the  target  plate  was  3  in  (7.62  cm),  one  half  of  the 
roof  outer  diameter.  The  structure  was  set  in  an  excavation  ensuring  that  its  center 
is  under  the  center  of  the  electromagnet  holding  the  impactor  ball.  Special  attention 
was  given  to  insure  leveling  of  the  roof  surface  and  sufficient  density  of  the  founding 
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sand  under  the  structure  floor.  Then  the  sand  backfill  was  poured  and  tapped  in 
layers  of  about  2  in  around  and  above  the  structure.  The  target  plate  was  placed  on 
the  soil  surface  and  centered  under  the  impactor  ball.  It  was  then  leveled  and 
tapped  to  provide  proper  contact  with  the  sand.  The  system  was  impacted  several 
times  to  assure  consistent  response  of  motion  and  pressure  records. 

22.4  Experiments  with  Reinforced  Micro-concrete  Roof  Slabs. 

2 2.4.1  Reinforced  micro-concrete  slabs. 

An  appropriate  way  to  model  reinforced  concrete  structural  elements  was 
recommended  in  research  by  Cunnigham,  Townsend,  and  Fagundo1121.  The 
preparation  and  treatment  of  the  mixture,  as  well  as  of  the  micro-reinforcement, 
were  performed  following  their  modeling  technique.  The  slabs  were  reinforced  with 
spot  welded  meshes.  The  mixture  for  the  micro-concrete  slabs  was  made  of  gypsum, 
fine  sand,  and  water.  The  use  of  gypsum  cement  provided  rapid  curing  time,  small 
particle  size,  and  low  distortion1371.  The  sand  was  sieved  into  three  portions,  as 
described  in  Figure  2-9.  The  average  weight  density  of  the  reinforced  slabs  was  130 
pcf  (20.4  kN/m3).  All  the  reinforcing  wires  were  deformed  to  obtain  similarity  to  the 
cold  deformed  conventional  reinforcement.  The  desired  deformation  was  obtained 
by  passing  the  wire  through  two  pairs  of  perpendicular  knurls1121.  The  wires  were 
preloaded  to  a  stress  near  the  yield  point  to  remove  any  initial  deformation.  After 
preloading  the  wires,  the  miniaturized  reinforcing  mesh  was  spot  welded.  Specimens 
of  the  wires  were  tested  in  tension  in  a  Material  Testing  System  (MTS)  machine  and 
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had  an  average  yield  strength  of  51000  psi  (352  MPa). 

The  compressive  strength  of  the  micro-concrete  mixture  was  tested  on  1  in  by  3 
in  (2.54  cm  by  7.62  cm)  cylinders,  and  the  micro-reinforcement  wires  were  tested  in 
simple  tension  tests.  These  tests  were  performed  on  an  MTS  on  samples  of  wires 
that  were  not  welded,  and  on  specimens  from  the  spot  welded  meshes,  to  verify  that 
the  spot  welding  technique  does  not  influence  the  properties  of  the  reinforcement. 
No  significant  difference  was  found  between  the  two,  in  their  yield  strength  or 
general  stress-strain  behavior.  The  properties  of  the  slabs  are  given  in  Table  2-1. 
The  forms  for  the  slabs  were  made  of  plexiglas.  Figures  2-10-a  through  2-10-d  show 
a  typical  process  of  casting  the  roof  slab  specimens. 

In  addition  to  the  strength  tests,  two  sets  of  structural  tests  were  performed  with 
the  circular  slab  specimens.  The  slabs  were  subjected  to  a  uniformly  distributed 
pressure  and  their  midspan  displacement  was  measured.  The  uniform  load  was 
applied  through  oil  pressure  in  a  steel  chamber  above  the  plate,  and  the  pressure  and 
displacement  were  measured  by  a  pressure  transducer  and  an  LVDT,  respectively 
(Figure  2-11).  In  the  first  set  of  tests  the  pressure  was  static,  and  the  tests  were 
controlled  by  a  constant  rate  deflection  of  the  slab  center.  The  purpose  of  the  static 
tests  was  to  ascertain  the  structural  behavior  of  the  small  scale  slabs  as  compared  to 
larger  similarly  reinforced  concrete  elements.  The  results  also  assisted  in  analyzing 
the  soil-structure  system.  The  slab  specimens,  similar  in  geometry  to  the  roof  slabs, 
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were  tested  under  controlled  conditions  and  were  observed  for  cracking, 
reinforcement  yielding,  and  ultimate  capacity.  This  study  enabled  a  better  analysis 
of  the  response  of  the  roof  slabs  (see  also  2.3.2  and  2.4.2).  The  experimental  system, 
and  the  results  of  the  static  tests  are  described  in  detail  by  Chen  et.  al.[10].  The 
second  set  of  experiments  were  performed  with  dynamic  loading  on  the  same  setup 
but  with  load  control  instead  of  displacement  control. 

It  was  found  in  the  static  tests  that  the  small  scale  specimens  behave  in  a 
comparable  manner  to  larger  scale,  similar  slabs.  The  impact  pressure  induced  by 
the  test  apparatus,  had  a  much  larger  duration  than  that  which  loaded  the  roof  in  the 
soil-structure  system.  It  was,  therefore,  not  possible  to  study  directly  the  uncoupled 
response  of  the  roof  slab  specimens  under  such  impact  loading.  Nevertheless,  the 
pressure  in  those  tests  had  a  rise  time  of  about  20  msec,  which  was  about  300  times 
faster  than  the  highest  rate  of  loading  in  the  static  tests.  Therefore,  the  mid-span 
displacement  record  was  used  in  the  evaluation  of  the  buried  roof  response,  and  as 
a  reference  in  the  numerical  analysis  (see  2.4.3  and  3.2.5). 

2.2A2  Boundary  conditions. 

With  the  assumption  that  the  plexiglas  walls  are  rigid  enough,  the  boundary 
conditions  of  the  roof  were  determined  according  to  its  means  of  support.  Two  types 
of  supports  were  used:  simply  supported  and  clamped.  The  simply  supported  fixture 
consisted  of  three  "L"  shaped  aluminum  clamps,  3/8  in  (0.95  cm)  wide,  which  were 
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attached  to  the  plexiglas  walls,  equally  spaced,  120  degrees  apart  from  each  other. 
The  top  portion  of  these  clamps  provided,  together  with  the  walls  of  the  structure, 
a  vertical  support  that  did  not  carry  any  moment  (figures  2-12-a  and  2-12-b).  The 
clamped  support  (figures  2-13-a  and  2-13-b)  consisted  of  an  aluminum  annular  cap 
which  was  mounted  over  the  slab  and  had  an  "L"  shaped  cross  section.  The  cap  had 
an  interior  diameter  identical  to  the  plexiglas  cylinder  (5  in)  and  was  attached  to  it 
by  6  steel  studs  that  prevented  any  vertical  movement  of  the  clamping  cap.  A  flat, 
1/32  in  (0.08  cm)  thick,  and  1/2  in  (1.27  cm)  wide,  steel  ring  was  placed  over  the 
roof  slab  and  under  the  clamping  aluminum  cap,  and  was  tightened  against  the  slab 
by  8  screws.  The  wall  and  top  part  of  the  clamping  cap  were  3/16  and  1/8  in  thick 
(0.48  and  0.32  cm),  respectively.  The  cap  was  rigid  enough  to  clamp  the  roof  (by 
preventing  a  rotation  of  its  edge)  but  not  too  large  or  thick  to  affect  the  general 
shape  of  the  structure. 

In  order  to  examine  the  effectiveness  of  the  supports  the  following  test  was 
conducted:  a  doubly  reinforced  slab  (#20  wire,  A, =9.5 115  -lO"4  in2  (0.0061  cm2),  @ 
1"  (2.54  cm)  in  each  direction),  0.43  in  (1.1  cm)  thick,  was  loaded  by  a  short, 
concentrated,  low  impact  at  its  center  (see  also  test  No.  6  in  Table  2-1).  The  free 
vibration  response  of  the  slab  was  recorded  to  obtain  its  fundamental  natural 
frequency.  If  the  clamped  boundary  fixture  is  effective,  it  is  expected  to  cause  a 
natural  frequency  higher  than  the  frequency  of  the  simply  supported  slab.  An  upper 
bound  for  the  natural  frequency  of  a  circular  slab  is  given  by  the  expression132,43491: 
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D_ 
P  h 


(2.11) 


where:  w  is  the  frequency  (rad/sec) 

D  is  the  stiffness  of  the  slab,  defined  as: 

D  =  E-h3  /  12(l-v2) 

E  is  Young’s  modulus 

h  is  the  thickness  of  the  slab  (0.43  in;  1.09  cm) 
v  is  Poisson  ratio 

p  is  the  mass  density  of  the  slab  (1.95 -104  lbs.-se^/in4;  2084  KgM/m4  , 
calculated  by  dividing  the  actual  weight  of  the  slab  by  its  volume  and  by  g) 
R  is  the  radius  of  the  slab  (2.5  in;  635  cm) 

a  is  a  constant  which  depends  on  the  boundary  conditions  and  on  the 
vibration  mode. 

For  the  case  of  a  fixed  circular  plate  a  =  10.22  for  the  lowest  mode  of  vibration. 
In  the  case  of  a  simply  supported  plate,  there  is  also  a  dependency  on  v  (e.g.,  Ref. 
32):  For  v  =  0.3,  and  0.25,  a  =  4.977,  and  4.858,  respectively.  Poisson’s  ratio  of  0.2 
was  assumed  for  the  micro-concrete. 


Figures  2-14-a  and  2-14-b  show  the  response  of  the  slab.  The  results  indicate  that 
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the  clamping  fixture  is  relatively  effective,  as  it  raised  the  natural  frequency  of  the 
slab  from  about  2800  rad./sec.  (446  Hz)  under  simple  support  conditions  to  about 
14000  rad./sec.  (2230  Hz)  when  it  was  clamped.  From  equation  2.11  and  the 
corresponding  values  of  a,  the  theoretical  values  of  a  simply  supported  slab  is 
expected  to  be  about  half  of  that  of  a  fixed  slab.  The  experimental  frequency  of  the 
slab  under  simply  supported  boundary  conditions  was  even  lower,  compared  to  that 
of  the  clamped  case.  However,  the  main  objective  of  this  part  of  the  experiment  was 
achieved  by  obtaining  two  different  types  of  boundary  conditions  which  generated 
two  different  natural  frequencies  for  the  embedded  roof.  It  should  be  noted  that 
while  the  theoretical  analysis  does  not  take  into  account  lateral  displacement  of  the 
slab  edge,  such  displacement  is  possible,  especially  when  a  concrete  slab  is  resting  on 
a  plexiglas  support.  Boundary  lateral  motion  decreases  the  stiffness  of  the  slab, 
thereby  decreasing  its  natural  frequency,  as  indeed  was  observed  in  the  above  test. 
Additional  indication  for  possible  lateral  motion  in  the  simply  supported  case  is  given 
in  the  numerical  analysis  (see  3.2.5). 

This  test  also  served  as  an  indirect  method  to  evaluate  the  initial  Young’s  modulus, 
E,  of  the  micro-concrete  with  the  above  parameters.  The  analytical  upper  bound  for 

the  lowest  natural  frequency  of  a  fixed  slab  is  15.18/E  rad/sec  •  The  measured 
frequency  value  of  14000  rad/sec  yields  a  relatively  low  Young’s  modulus  of  0.89106 
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psi  (6137  MPa).  Similar  measurements  were  taken  for  the  slabs  of  tests  7  and  8  and 
yielded  corresponding  values  of  1.05 -106  and  0.98 -106  psi  (7240  and  6757  MPa), 
respectively.  Values  of  Young’s  modulus  which  were  measured  in  compression  tests 
of  micro-concrete  indicated  higher  values  but  also  showed  considerable  scattering. 
It  is  interesting  to  note  that  a  similar  small  value  of  elastic  modulus  for  micro¬ 
concrete  was  also  reported  by  Shin[40]. 

22.43  Test  procedure. 

The  configuration  of  the  test  setup  was  similar  to  that  of  the  experiment  with  the 
plexiglas  roof  (see  section  2.2.3).  In  tests  1  to  5a  the  roof  was  embedded  3  in  (7.62 
cm)  under  the  soil  surface  and  in  the  other  tests  the  DOB  was  decreased  to  one  half 
of  the  roof  clear  span,  i.e.  2.5  in  (6.35  cm).  The  boundary  conditions  of  the  roof 
were  described  in  the  previous  section  and  are  given  for  each  test  in  Table  2-1.  In 
the  tests  where  the  roof  slab  was  simply  supported,  the  load  was  measured  at  the 
center,  2  in  (5.08  cm)  from  the  center,  and  2.5  in  (6.35  cm)  from  the  center.  This 
configuration  of  load  measurement  was  not  possible  in  the  case  of  the  clamped  slabs, 
because  the  inner  diameter  of  the  clamping  cap  was  5  in  (12.7  cm)  and  each  load  cell 
had  a  diameter  of  0.5  in  (1.27  cm).  Therefore,  in  tests  3  through  5b  and  in  test  8  the 
load  cells  were  located  at  the  center  and  2  in  (5.08  cm)  from  the  center. 

The  steel  ball  impactor  was  dropped  from  an  initial  height  of  25  in  (63.5  cm)  or 
50  inches  (127  cm),  and  the  center  displacement  signal  was  examined  after  each  drop 
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for  residual  plastic  deflection,  indicating  possible  failure.  When  such  indication  was 
observed,  the  structure  was  excavated  and  examined.  Table  2-2  lists  the  sequence 
of  hits  in  the  tests.  As  a  consequence  of  the  preliminary  test  with  the  plexiglas  roof, 
the  target  plate  was  pressed  and  tapped  into  the  soil  surface  to  avoid  a  relatively 
large  settlement  under  the  first  impact. 
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23  Results. 

2.3.1  Experiment  with  an  Elastic  Plexiglas  Roof  Slab. 

It  was  observed,  in  both  free  field  experiments,  and  the  experiments  with  the 
structure,  that  under  the  first  impact  there  was  a  relatively  significant  settlement,  up 
to  0.5  inch  (1.27  cm),  of  the  target  plate.  Therefore,  the  initial  cover  was  0.5  inch 
deeper  than  the  test  design  cover,  and  the  results  of  the  first  hit  were  usually 
disregarded.  This  approach  was  confirmed  to  be  appropriate,  by  measuring  a 
constant  depth  of  cover,  and  similar  pressure  and  motion  records  in  the  results 
succeeding  the  second  hit,  for  each  drop  of  the  impactor.  It  was  also  reasonable  to 
anticipate  this  behavior,  since  the  plexiglas  slab  neither  failed  nor  yielded  in  any  of 
the  tests  performed  at  this  stage. 

Figure  2-15  shows  the  displacement  records,  w,,  w2,  and  w3.  Theoretically,  a 
double  integration  of  the  acceleration  records  should  provide  the  absolute  vertical 
displacement  of  the  center  of  the  roof  (denoted  as  w2)  and  of  the  whole  (rigid) 
structure  (denoted  as  w3),  such  that  Wj  =w2-w3,  where  Wj  is  the  LVDT  output  (see 
figure  2-8).  As  shown  in  Appendix  B  the  LVDT  signal  is  more  accurate  than  that 
of  the  doubly  integrated  accelerometer  signal  but  slower  to  respond.  Since  the 
records  taken  during  the  first  100  nsec  of  loading  are  important,  the  signal  of  the 
accelerometer  serves  to  indicate  the  beginning  of  the  motion  of  the  structure. 
Accordingly,  Table  2-3  describes  the  data  acquisition  method  for  these  experiments. 
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Figure  2-16  shows  the  pressure  and  corresponding  velocity  signals  as  were  recorded 
in  the  test,  over  a  period  of  4  msec.  As  can  be  seen,  the  significant  portion  of  these 
signals  occurred  only  within  the  first  2  msec  of  the  test,  where  the  time  of  impact 
initiation  under  the  target  plate  center  was  set  to  0.0  sec.  Therefore,  the  response 
over  that  period  was  magnified  and  shown  in  figures  2-17  and  2-18.  The  vertical 
stresses,  together  with  the  velocity  and  displacement  of  the  roof,  are  shown  in  figure 
2-17.  Figure  2-18  demonstrates  the  difference  between  the  stress  distribution  in  the 
free  field  and  over  the  buried  structure.  The  impulse  records,  integrated  from  the 
pressure  signals,  are  given  in  figure  2-19. 

2 32  Experiments  with  Reinforced  Micro-concrete  Roof  Slabs. 

The  roof  slabs  which  were  buried  under  3  inches  (7.62  cm)  of  sand  (tests  1  through 
5,  in  Tables  2-1  and  2-2)  showed  slight  or  no  visible  damage.  The  most  visible 
damage  was  observed  in  tests  1  and  2.  Under  an  external  impact,  which  was 
generated  by  a  height  of  drop  of  50  inches  (127  cm),  the  slab  was  cracked  and  the 
LVDT  was  detached.  The  cracks  were  at  the  bottom  of  the  plate,  and  had  a  radial 
shape.  These  cracks  were  the  only  damage  caused  and  indicate  that  the  response 
was  flexural.  In  tests  3,  4  and  5a  there  was  only  slight  radial  cracking,  even  under 
external  impact  from  a  height  of  drop  of  90  inches  (228.6  cm). 

In  several  cases  the  signal  of  the  LVDT  indicated  relatively  large  displacements. 
When  the  structure  was  excavated  no  visible  damage  was  observed,  but  the  core  of 
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the  LVDT  was  detached  from  the  roof  slab.  Due  to  the  high  acceleration  of  the  roof 
center  (up  to  5000  g’s)  the  fixture  which  held  the  core  and  was  glued  to  the  slab, 
became  detached,  although  it  had  a  relatively  low  mass  (0.0178/g  lb/g;  8.1  gr).  The 
detachment  of  the  LVDT  fixture  also  indicates  transmission  of  a  tension  wave  from 
the  bottom  surface  of  the  roof,  (see  also  the  discussion  of  the  results  in  the  next 
section). 

Similar  types  of  response  and  damage  to  the  roof  slabs  were  observed  in  tests  5b 
through  8,  where  the  DOB  was  decreased  to  2.5  in  (6.35  cm).  A  permanent  residual 
deflection  of  the  roof  center  was  observed  after  the  impactor  ball  was  dropped  from 
a  height  of  90  in  (228.6  cm),  except  for  test  7,  where  initial  cracks 
appeared  under  an  impact  generated  by  a  height  of  drop  of  25  in  (63.5  cm),  more 
cracking  was  seen  under  an  impact  from  a  drop  of  75  in  (190.5  cm). 

Structural  damage  to  the  slabs  was  indicated  by  the  appearance  of  visible  cracking 
and  plastic  residual  displacement.  Slight  to  mild  degree  of  damage  of  the  roof  slab 
was  observed  in  all  tests,  but  there  was  no  catastrophic  failure.  In  general,  the 
clamped  slabs  showed  higher  stiffness,  which  was  demonstrated  by  higher  frequency 
response,  and  relatively  lower  amplitude  of  the  displacement  and  velocity. 
Nevertheless,  no  cracks  were  detected  on  top  of  the  clamped  slabs,  which  indicates 
that  these  slabs  were  not  fully  fixed  at  their  supports,  section  2.2.4.2.  The  maximum 
measured  central  deflection  in  the  tests  was  only  about  0.07  in  (0.18  cm)  in  test  6, 


under  the  highest  drop  of  the  impactor  ball  (see  also  results  of  test  6,  4th  hit,  in 
Appendix  D).  These  relatively  low  displacements  also  explain  the  mild  level  of 
damage  which  was  observed  (see  also  Table  2-2  and  the  discussion  in  section  2.4.3). 

Table  2-2  lists  the  sequence  of  the  hits  in  the  tests  that  were  performed  and  their 
intermediate  and  final  observations.  A  typical  result,  which  includes  measurements 
of  the  pressure  at  the  roof-soil  interface,  and  of  the  roof  motion,  are  given  in  figures 
2-20  through  2-23  and  discussed  in  the  following  sections.  Other  results  are  shown 
in  Appendix  D. 
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2A  Discussion. 

2.4.1  Physical  Mechanisms  and  Phenomena  that  are  Expected  to  be  Involved. 

Analysis  and  discussion  of  the  test  results  are  given  in  view  of  physical  mechanisms 
and  phenomena  that  are  expected  to  be  involved  in  the  system  and  were  discussed 
or  reported  in  the  literature  (see  also  Chapter  1).  They  are: 

1.  Incident  free  field  stress  distribution,  which  is  the  stress  field  in  the  soil  under 
external  load  with  no  embedded  structure.  The  free  field  pressure  constructs 
the  initial  loading  on  the  embedded  structure,  whose  response  causes  deviation 
of  the  interface  load  from  the  free  field  pressure  at  the  same  location  and  time. 

2.  The  initial  total  pressure  acting  on  the  structure  is  the  incident  pressure,  p'  plus 
the  reflected  pressure  from  the  soil-structure  boundary,  pr.  It  is  determined  by 
the  boundary  conditions  at  the  soil-structure  interface  and  the  effect  of  the 
mechanical  impedance  of  the  soil.  The  mechanical  impedance  is  a  measure  of 
the  resistance  of  a  system  or  medium  to  an  applied  load.  In  a  single  degree  of 
freedom  system  the  impedance  is  equal  to  pc,  where  p  is  the  density  and  c  is 
the  longitudinal  wave  velocity.  As  long  as  there  is  no  gap  between  the  soil  and 
the  structure,  and  for  a  plane  pressure  wave,  the  impedance  gives  the  relation 
between  the  pressure,  pf  and  the  velocity  of  the  soil  particles,  v^,: 
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pl  =  pf  =  (pc)^  (2.12) 

The  reflected  pressure  wave  and  the  total  pressure,  p,  at  the  soil-structure 
interface  are  given  by: 

pT  =  (pcXv^  -  vj  =  pf  -  (pc)v^  (2.13) 

P  -  Pl  +  Pr  =  2//  -  (pc)^  (2.14) 

where  vstr  is  the  velocity  of  the  structure.  (vstr  is  formally  the  velocity  of  the 
structure  particles  at  the  interface.  For  a  relatively  rigid  roof,  such  as  concrete, 
and  an  incident  wave  with  a  relatively  long  duration,  such  as  the  free  field 
impulse,  vstr  is  the  actual  velocity  of  the  roof). 

It  is,  therefore  noted  that  this  mechanism,  associated  with  wave  propagation 
phenomena,  is  dominated  by  the  boundary  velocity  of  the  particles  of  the 
medium  in  which  the  wave  propagates.  When  the  medium  is  bounded  by  an 
elastic  boundary  its  particle  velocity  is  also  the  boundary  velocity,  unless  there 
is  a  gap.  In  the  present  system,  the  structure  constitutes  a  boundary  for  the  soil 
particles  which  are  resting  on  it  until  motion  starts.  When  loaded  by  an  impact, 
the  roof  responds  by  gaining  velocity,  which  also  affects  the  interface  pressure. 
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Equation  2.14  implies  that  in  the  initial  time  of  loading,  when  the  velocity  of  the 
structure  is  not  yet  developed,  the  incident  pressure  is  amplified  at  the  soil- 
structure  interface.  When  later  structure  velocity  is  developed  in  response  to  the 
applied  load,  the  soil-interface  pressure  is  affected  by  the  interface  velocity- 
pressure  relation  according  to  equation  2.14  and  is  expected  to  be  reduced. 

3.  Tension  waves  normal  to  the  roof  plane  are  generated  when  the  portion  of  the 
propagating  compression  impact,  transmitted  into  the  slab,  arrives  at  the 
bottom  face  of  the  roof. 

4.  Arching  in  the  soil  is  expected  with  developing  structural  displacement.  By 
transferring  the  pressure  from  the  more  flexible  parts  of  the  structure  to  its 
stiffer  parts,  the  soil  interacts  with  the  deflecting  structure,  mainly  through  its 
shear  capacity.  This  mechanism  is  dominated  by  displacements  of  the  structure 
and  of  the  soil.  In  case  of  relatively  large  displacement  of  the  roof  center 
arching  tends  to  decrease  the  pressure  over  this  region  and  increase  it  over  the 
edge.  In  case  of  rigid  body  settlement  of  the  structure  the  DOB/span  ratio 
decreases.  Consequently  the  arching  action  over  the  roof  is  expected  to 
decrease  and  the  pressure  at  the  center  of  the  roof  may  increase  again.  The 
phenomenon  of  arching  is  also  explained  and  described  in  Chapter  1. 


5.  Possible  soil-structure  gap  over  the  center  ("Double  hit").  With  sufficient 
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structure  velocity  a  gap  might  develop  between  the  soil  particles  and  the 
moving  structure.  Such  gap  would  be  indicated  by  decrease  of  the  interface 
pressure  to  a  minimum  of  zero.  When  the  center  velocity  decreases,  the 
velocity  of  the  sand  above  this  region  may  again  match  the  velocity  of  the  roof 
and  reload  it. 

2.4 2  Sequence  of  the  Phenomena  in  the  Response  of  the  Soil-Structure  System. 

The  results  are  examined  in  light  of  the  physical  mechanisms  listed  in  the  previous 
section  and  phenomena  which  are  expected  to  be  involved  in  the  system  response. 
The  following  observations  are  of  the  results  of  test  number  8,  5th  hit,  which  was 
performed  with  a  clamped  slab,  and  a  DOB  of  2.5  inches  (Table  2-1).  Results  of 
other  tests,  given  in  Appendix  D,  showed  similar  observations. 

Comparison  of  the  pressure  records  over  the  roof  center  and  edge  to  the 
corresponding  free  field  records  shows  that  after  arriving  at  the  roof  level  the 
pressure  initially  follows  the  free  field  behavior.  The  initial  free  field  loading  of  the 
roof  was  indeed  measured  and  recorded  in  all  the  tests  and  lasted  for  0.02  to  0.05 
msec  (0.15<t<0.20  msec  in  figure  2-20).  This  period  of  time  was  observed  to  be 
insensitive  neither  to  the  two  DOB  which  were  used  nor  to  the  boundary  conditions; 
it  increased  in  the  test  with  the  plexiglas  roof  (about  0.1  msec  from  the  beginning  of 
loading  in  figure  2-18).  The  rise  time  to  the  peak  was  also  relatively  longer  under 
the  highest  external  impact  in  test  6,  and  in  test  8.  Accordingly,  the  displacements 
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under  the  longer  and  higher  center  impacts  were  larger  than  in  the  other  hits  (e.g., 
deflection  of  0.07  in  (0.18  cm)  in  the  4th  hit  in  test  6  versus  maximum  deflections  of 
about  0.03  in  (0.08  cm)  under  the  other  hits  of  that  test,  Appendix  D). 

No  amplification  of  the  interface  pressure,  relative  to  the  free  field  level  was 
observed.  The  velocity  which  is  required  to  prevent  an  amplification  of  the  incident 
pressure  can  be  evaluated  from  setting  the  interface  pressure,  p,  equal  to  the  free 
field  pressure,  pf  in  equation  2.14.  The  relevant  properties  of  the  soil  are: 

Psou  =  161 0^  U>-secz/in4  (1710  KgM/m3) 

cu  =  - — —  =  15000 in/sec  (381  m/sec) 

arrival  time 

For  the  condition  p=pf  to  be  satisfied  at  a  level  of  100  psi  (0.69  MPa),  the  roof 
velocity  has  to  be  in  the  order  of  40  in/sec  (1.02  m/sec).  However,  during  the  initial 
loading  period  the  velocity  of  the  structure  was  developing  from  zero  (rest)  to  less 
than  50  in/sec  (1.27  m/sec;  0.15  msec<t<0.2  msec  in  figure  2-20). 
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In  the  apparent  lack  of  sufficient  roof  velocity  the  interface  pressure  readings  may 
be  explained  by  considering  a  tension  wave  reflected  from  the  bottom  of  the  roof. 
The  velocity  of  an  elastic  wave  through  the  roof  thickness  is  given  by  equation  2.15: 
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(2.15) 

where:  cL  is  the  longitudinal,  or  compressive  wave  velocity. 

E  is  the  concrete  Young’s  Modulus, 
v  is  the  concrete  Poisson  ratio, 
p  is  concrete  mass  density. 

The  mass  density  of  micro-concrete  was  1.95  10"4  lb-sec^/in4  (2084  KgM/m4). 
Poisson’s  ratio  of  concrete  can  be  assumed  as  0.1  to  0.2,  which  in  either  case  does 
not  influence  significantly  the  value  of  the  velocity  given  by  equation  2.15.  Table  2-4 
shows  the  wave  speed  values  for  these  material  properties,  for  Young’s  modulus  of 
2-106  psi  (13790  MPa)  and  5-106  psi  (34475  MPa),  and  the  corresponding  time,  tT 
it  would  take  for  a  reflected  tension  wave  to  arrive  at  the  upper  face  of  the  roof. 
The  time  tf  is  calculated  according  to  equation  2.16,  assuming  that  the  velocity  of  the 
incident  compression  wave  and  of  the  reflected  tension  wave  are  the  same: 

.  =  iJL  (2.16) 

'l 

where  h  is  the  thickness  of  the  roof. 


46 

As  can  be  seen  in  Table  2-4,  from  the  incidence  of  the  compression  impact  on  the 
roof  upper  face,  up  to  0.01  msec  are  required  for  a  tension  wave  to  be  reflected  and 
arrive  from  the  bottom  free  face  of  the  roof.  The  test  results  show  that  until  that 
time  the  free  field  pressure  was  still  low  enough  and  the  measurements  of  the 
interface  pressure  were  likely  to  have  a  relatively  high  level  of  noise  to  signal  ratio. 
Thus,  it  is  possible  that  even  if  it  existed,  an  amplification  may  have  not  been 
detected. 

A  reflected  tension  wave  might  have  caused  a  decrease  of  the  soil-roof  pressure 
when  propagating  back  to  the  upper  roof  surface,  and  spalling  at  the  bottom  face  of 
the  roof.  No  spalling  was  observed  in  any  of  the  tests,  which  is  reasonable  since  the 
peak  pressure  value,  measured  under  the  highest  drop  of  the  impactor,  was  about  160 
psi  (1.1  MPa;  figure  2-20).  An  equal  tension  stress  wave  reflected  from  the  free 
bottom  of  the  roof  could  be  sustained  by  the  concrete.  An  observation  which  might 
resemble  spalling  and  indicate  a  reflected  tension  wave  at  the  roof  bottom,  was  the 
local  detachment  of  the  LVDT  fixture  at  the  slab  center,  which  occurred  several 
times  (see  also  Table  2-2). 

About  0.02  to  0.05  msec  after  hitting  the  roof  both  the  pressure  over  the  center 
and  over  the  edge  started  to  drop.  In  comparison  with  corresponding  free  field 
reference  signals  this  drop  was  a  deviation  caused  by  the  roof  gaining  velocity  upon 
being  hit.  Therefore,  it  was  concluded  that  for  the  first  0.05  to  0.1  msec  the  response 
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of  the  roof  involved  a  mechanism  appropriate  to  wave  propagation  in  a  medium  with 
an  elastic  boundary,  where  the  velocity  of  the  boundary  affects  the  normal  pressure 
acting  on  it.  The  time  of  initial  deviation  from  the  free  field  pressure  is  marked  with 
"A"  in  the  figures. 

Approximately  0.1  to  0.15  msec  from  the  initial  loading  of  the  roof,  the  pressure 
over  the  center  continued  to  drop  while  over  the  roof  edge  it  rose  again  and  returned 
to  the  free  field  level  (Marked  "B"  in  figure  2-20,  at  t*0.24  msec),  although  the  edge 
was  still  moving  downwards  with  a  negative  velocity  (figure  2-21,  t«0.24  msec).  This 
part  of  the  response  indicates  the  development  of  arching  identified  by  a  decrease 
of  the  normal  pressure  over  the  more  flexible  central  part  of  the  roof  and  its  increase 
over  the  stiffer  edge,  due  to  relative  displacement  of  the  center.  Note  that  in  the 
axisymmetric  geometry  of  this  system,  the  arching  effect  is  more  likely  to  occur  as  a 
dome  action  in  the  soil.  It  is  still  referred  to  as  arching  because  it  is  the  usual 
terminology  for  this  mechanism. 

The  pressure  over  the  center  dropped  to  zero  (or  maintained  a  moderately  higher 
minimum  in  other  tests),  while  the  pressure  record  over  the  edge  of  the  roof  did  not 
drop  to  zero  during  the  first  1  to  1.5  msec  of  the  response.  This  phenomenon 
indicates  that  in  some  cases  where  the  pressure  over  the  center  dropped  to  zero,  the 
velocity  of  the  roof  was  sufficient  to  cause  detachment  of  its  center  part  from  the 
sand  above  it  for  a  period  of  up  to  0.2  msec,  as  also  can  be  seen  in  the  impulse 
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records  (e.g.,  points  "B"  to  MC"  in  figures  2-20,  2-21,  2-22),  which  show  a  decrease  in 
the  slope  of  the  center  impulse.  The  period  of  minimum  or  zero  slope  indicate 
minimum  pressure,  which  when  equal  to  zero  can  be  interpreted  as  detachment  of 
the  soil  particles  from  the  roof.  The  minimum  pressure  over  the  roof  center  is 
marked  "C"  in  the  figures. 

About  0.17  to  0.20  msec  after  the  initial  loading  of  the  roof  (t»0.35  msec  in  figures 
2-20  and  2-21),  and  during  a  period  of  about  0.20  msec,  there  was  a  second  increase 
of  the  normal  pressure  over  the  center  of  the  roof,  while  the  pressure  over  the  edge 
slightly  decreased  below  the  level  of  the  free  field.  During  this  period  of  time  the 
motion  of  the  roof  center  was  still  downward  (negative  center  velocities  during  t  < 
0.58  msec  in  figure  2-21).  The  amplitude  of  the  center  velocity  decreased  at  that 
time,  and  the  velocity  of  the  sand  above  this  region  may  have  again  matched  the 
velocity  of  the  roof  and  reloaded  it.  It  was  also  possible  that  due  to  a  rigid  body 
displacement  of  the  structure,  causing  the  displacement  of  the  edge  of  the  roof, 
negative  arching  (or  a  decrease  in  the  previously  developed  arching)  was  developed 
over  it  and  the  pressure  over  the  center  increased. 

Further  analysis  shows  that  0.3  to  0.4  msec  after  the  arrival  of  the  incident  impact, 
the  velocity  of  the  roof  center  became  positive  (Figure  2-21  at  t  =  0.58  msec,  marked 
"D"  in  the  figures),  indicating  an  upward  motion,  and  the  normal  pressure  over  that 
part  of  the  roof  increased  again  in  a  similar,  but  opposite,  manner  to  the  decrease 
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of  the  pressure  at  the  beginning  of  the  response  of  the  structure,  at  t  «  0.2  msec. 
Here,  however,  the  increase  in  the  pressure  was  not  observed  in  all  the  cases  (e.g., 
see  results  of  test  2,  4th  hit  in  Appendix  D). 

2.43  Influence  of  the  Involved  Phenomena  on  the  Soil-Structure  Response. 

The  physical  mechanisms  indicated  by  the  experimental  observations,  occurred  at 
the  sequence  described  in  the  previous  section.  Their  influence  on  the  response  of 
the  structure  is  evaluated  by  examination  of  the  roof  displacement.  Inspection  of  the 
cracks  of  the  roof  slab  of  test  8  after  the  fifth  hit  and  a  recorded  residual 
displacement  of  about  0.02  in  (0.05  cm)  of  its  center,  indicated  a  state  of  flexural 
damage. 

More  information  about  the  behavior  of  the  slab  is  obtained  from  examination  of 
typical  load-deflection  curves  of  circular  micro-concrete  slabs  of  similar  thickness  and 
span.  Tests  with  uniform  pressure  were  performed  on  slabs  that  were  fixed  to 
rotation  and  to  lateral  displacement  at  the  boundary1101.  The  pressure  was  applied 
with  center  displacement  control  at  a  slow  rate  of  0.0133  in/min  (0.034  cm/min)  and 
at  faster  rate  of  0.52  in/min  (1.32  cm/min),  and  had  a  rise  time  to  peak  pressure  of 
about  900  sec  and  6  sec,  respectively.  Additional  tests  were  performed  under 
pressure  control,  that  had  a  high  rise  time  of  about  20  msec.  The  slab  was  simply 
supported  to  allow  later  convenient  numerical  modeling  (see  also  Chapter  3).  Figure 
2-23  shows  the  load-deflection  curves  of  those  tests.  The  rise  time  of  the  first  central 
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pressure  peak  in  test  8  (5th  hit)  was  0.05  msec.  The  pressure-time  records  (over  the 
center  and  over  the  edge)  were  plotted  with  respect  to  the  displacement-time  record 
eliminating  the  common  time  parameter,  to  give  the  load-displacement  history  of  that 
test  (figure  2-23).  As  indicated  by  figure  2-23,  the  micro-concrete  slabs  in  the  load- 
displacement  tests  sustained  a  maximum  load,  which  increased  at  decreasing 
displacement  with  increasing  rate  of  loading.  They  also  indicate  that  the  central 
impact  on  the  slab  which  lasted  less  than  0.1  msec  (figures  2-20  and  2-21)  would  have 
caused  more  damage  with  a  longer  duration,  such  as  that  of  the  free  field  impact 
(figure  2-20). 

This  point  is  also  demonstrated  in  figure  2-22  where  the  pressure  records  over  the 
structure  were  integrated  to  give  the  impulses.  Points  "A"  to  "D",  defined  in  the 
previous  section,  are  marked  in  the  figure.  They  show  the  initial  increase  of  the  load 
up  to  point  "AH,  where  the  structure  velocity  followed  by  displacement  (point  "B") 
caused  decrease  of  the  load  over  the  center  roof  up  to  point  "C",  During  this  time 
("A"  to  "C")  the  slope  of  the  impulse  over  the  center  dropped  to  a  plateau  which 
denotes  reduction  of  the  load. 

Therefore,  it  is  concluded  that  the  amplification  of  the  incident  pressure  wave  due 
to  the  relative  stiffness  of  the  structure  was  prevented,  by  the  tension  wave  reflected 
from  the  bottom  of  the  slab.  A  downward  velocity  and  displacement  of  the  roof 
further  decreased  the  load  acting  on  it.  Due  to  the  reduction  of  the  load  over  the 


51 


central  portion  of  the  roof  limited  damage  was  suffered,  and  further  catastrophic 
failure  was  prevented.  Consequently,  an  important  consideration  would  be  to 
evaluate  the  capacity  of  the  buried  structure  for  dynamic  impact  loading.  If  that 
capacity  is  higher  than  the  initial  portion  of  the  free  field  load,  then  the  other 
mechanisms  described  above  are  likely  to  decrease  the  free  field  load  so  that  the 
structure,  which  might  have  failed  under  a  load  equal  to  the  total  free  field  load,  is 
left  intact. 

Hence,  it  seems  that  the  portion  of  the  free  field  impulse,  acting  on  the  structure 
until  it  decreases  by  the  mechanisms  described  above,  has  an  important  effect  on  the 
response  or  even  capacity  of  the  structure  and  is  a  valid  result,  provided  that  the 
tension  wave  reflected  from  the  free  bottom  of  the  slab  is  either  low  enough  or  can 
otherwise  be  sustained  by  the  slab. 

2.4.4  Influence  of  the  Boundary  Conditions. 

The  effect  of  the  boundary  conditions  is  illustrated  by  comparing  the  response  of 
the  clamped  roof  in  test  3  with  that  of  the  simply  supported  roof  in  test  2.  Figure 
2-24  shows  the  external  pressure  under  the  center  of  the  target  plate,  as  a  result  of 
dropping  the  impactor  from  a  height  of  25  in  (63.5  cm)  in  tests  2  and  3.  The  slabs 
in  these  tests  had  a  similar  thickness  and  were  embedded  3  in  (7.62  cm)  under  the 
target  plate.  Under  a  similar  external  impact  (Figure  2-24)  the  clamped  slab 
exhibited  lower  maximum  displacement  (Figure  2-25;  the  displacement  curve  of  test 
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3  was  shifted  to  allow  convenient  comparison). 

When  beginning  to  analyze  the  soil-structure  interaction  aspects,  one  notes  the 
displacement  of  the  roof  slab,  which  effects  the  development  and  amount  of  arching, 
took  place  about  0.1  msec  after  the  beginning  of  the  structural  response.  Until  that 
time  the  load  was  not  affected  by  the  structure  boundary  conditions.  The  effect  of 
the  boundary  conditions  was  observed  when  the  response  of  the  structure  became 
significant.  The  impulse  records  over  the  roof  edge  and  center  in  the  two  cases  are 
given  in  figure  2-26  and  shows  lower  impulse  over  the  roof  center  in  the  simply 
supported  case,  and  a  corresponding  higher  impulse  over  the  edge  (dashed  and  solid 
lines  for  the  simply  supported  and  clamped  cases,  respectively,  in  figure  2-26).  This 
response  indicates  stronger  arching  over  the  more  flexible  simply  supported  roof.  It 
should  also  be  noted  that  although  not  measured,  it  is  reasonable  to  assume  a  more 
concave  deflection  shape  of  the  simply  supported  roof,  which  would  tend  to  increase 
the  more  flexible  zone  of  the  roof  and  decrease  the  load  over  it,  by  arching. 

It  seems  that  on  repeated  hits  the  stiffness  of  the  clamped  slabs  decreased,  either 
by  micro-cracking  or  by  some  loosening  of  the  grip  of  the  clamping  cap.  Noting  that 
the  decrease  of  the  pressure  over  the  roof  center  occurred  in  the  first  0.1  to  0.2  msec 
of  the  structural  response,  it  can  be  seen  that  the  roof  did  not  reach  its  maximum 
deflection  at  that  time.  As  a  result,  the  displacement  of  the  roof  center  in  the  first 
0.1  msec  of  the  response  was  larger  after  it  was  hit  repeatedly  (figure  2-28)  and 
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consequently  the  central  part  of  the  slab  was  less  loaded  than  in  the  previous  hits 


(compare  first  and  third  hits  in  figure  2-29). 

2.4.5  Effect  of  DOB. 

The  free  field  pressure  under  the  same  surface  impact  is  expected  to  undergo  more 
attenuation  at  deeper  depth.  At  shallow  depths  of  burial  the  relative  effect  of  other 
soil-structure  mechanisms  become  more  significant.  The  results  indicate  stronger 
arching  level  at  deeper  DOB.  Comparing  test  2  and  test  3,  which  were  done  at  DOB 
of  3  in  (7.62  cm)  and  2.5  in  (6.35  cm),  it  can  be  seen  in  figure  2-30,  that  at  the  time 
when  arching  occurred  (after  the  center  impulse  reached  a  plateau)  the  edge  to 
center  impulse  ratio  was  higher  in  test  2  where  the  DOB  was  deeper: 


^ edge 

_  0.008 

^edge 

_  0.011 

I center 

(■0.43msec  0.004 

I center 

(■0.42m*c  0.008 

DOB =3" 

DOB =2-5" 
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2t£  The  Value  of  the  Results  from  the  Small  Scale  Experimental  System. 

The  results  of  the  experiments  that  were  performed  with  the  small  scale  system 
have  two  main  applications: 

1.  Obtaining  detailed  information  that  will  enable  close  study  of  the  basic 
mechanisms  which  are  involved  in  a  soil-embedded  structure  system. 

2.  Enable  predictions  of  the  response  of  similar  larger  scale  systems. 

As  for  the  first  application,  it  is  sufficient  to  analyze  the  parameters  and 
experimental  records  as  they  are.  This  study  provides  an  analysis  of  the  small  scale 
system.  Nevertheless  it  is  important  to  evaluate  the  value  of  the  results  obtained 
here  for  the  understanding  and  analysis  of  larger  scale  systems  of  structures 
embedded  in  soil. 

For  this  reason  ,  dimensional  analysis  was  performed  based  on  the  theory  of 
similitude,  which  is  detailed  in  Appendix  C. 

Table  2-5  shows  the  ratios  that  should  be  used  when  extrapolating  any  variable 
from  the  small  scale  model  system  to  a  larger  prototype  system.  These  extrapolation 
rules  are  based  on  the  following  assumptions: 
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1.  The  geometrical  scaling  factor  is  equal  to  a  known  chosen  factor,  n. 


2.  The  materials  which  are  used  in  the  model  are  similar  in  their  mechanical  and 
weight  properties  to  those  which  are  used  in  the  prototype. 


A  well  known  conflict  is  introduced  when  similar  materials  are  used  for  the 
prototype  and  for  the  model,  and  the  gravitational  acceleration  is  not  scaled  (e.g.  by 
using  a  centrifuge)115,231.  This  causes  some  distortion  in  the  model  to  prototype 
extrapolation.  That  distortion,  however,  is  relatively  negligible  and  falls  within  the 
range  of  other  experimental  errors,  when  the  gravitational  effect  is  minor.  As 
demonstrated  in  Appendix  C,  if  a  variable  in  the  system,  "VAR",  is  assumed  to 
consist  of  a  static  portion,  "VAR*,",  and  of  a  dynamic  portion,  then  the  distortion 
depends  on  the  geometrical  scaling  factor  and  on  the  ratio  of  the  total  amplitude  of 
the  variable  to  its  static  portion,  according  to  equation  2.23: 


%  Distortion 


VAR1™- VAR. m  _  100 

VAR1™  ,  l  VARm 

l\n-l)VARmM 


(2.23) 


where:  VARtn>e  is  the  true  prototype  value 

VAR  is  the  extrapolated  prototype  value 

VARm  is  the  model  value  (usually  measured  in  the  test) 
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VAR,^,  is  the  static  portion  of  the  model  value 

In  the  case  of  a  dominant  dynamic  portion  the  distortion  decreases  as  it  does  when 
the  scaling  factor  is  low.  It  is  also  shown  in  Appendix  C  that  the  levels  which  were 
obtained  of  the  pressure  over  the  roof,  and  of  the  mid-roof  displacement  were  larger 
than  their  corresponding  static  values.  Hence,  a  possible  distortion  due  to  the 
unsealed  gravitation  is  low  to  moderate,  and  the  ratios  given  in  Table  2-5  are 
sufficient  for  a  model  to  prototype  predictions  based  on  the  current  system  with  a 
geometrical  scaling  factor  as  high  as  60. 

It  should  be  noted,  however,  that  the  external  pressure  cannot  be  extrapolated,  and 
that  the  relatively  low  level  impact  which  was  produced  in  the  current  system  (up  to 
500  psi;  3.45  MPa)  might  not  be  typical  of  a  full  scale  explosive  impact.  However, 
the  time  is  scaled  up  by  a  factor  of  n,  expanding  the  extrapolated  prototype  time  to 
a  longer,  more  realistic  duration. 

There  are  three  more  aspects  which  present  possible  difficulties  in  the 
extrapolation  of  the  small  scale  system  response: 

1.  The  response  of  the  reinforced  micro-concrete  in  the  model:  The  usage  of  this 
particular  technique  followed  the  recommendations  of  Cunningham,  Townsend, 
and  Fagundo'125  who  found  it  suitable  for  modeling  larger  scale  reinforced 
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concrete  elements.  Further  confirmation  was  obtained  in  the  observations  of 
Chen  et.  al.[10]  who  tested  similar  small  scale  slabs  and  found  their  response  to 
be  similar  to  the  response  of  larger  scale  slabs. 

2.  Modeling  of  the  backfill:  On  one  hand  the  use  of  20-30  Ottawa  sand  cannot 
model  geometrically  a  typical  soil  backfill  which  is  used  over  larger  scale 
structures.  On  the  other  hand  its  main  properties,  i.e.  zero  tensile  strength, 
shear  capacity,  and  granular  structure,  are  assumed  to  be  sufficient  for 
representing  the  primary  characteristics  of  such  backfill. 

3.  Strain  rate  effect,  and  size  effect15,61  might  be  significant  in  the  analysis  of  the 
response  of  concrete  elements.  This  is  especially  important  in  failure  initiating 
crack  propagation  in  quasi-brittle  materials  such  as  concrete111,1’ J.  For  the 
reinforced  concrete  slabs  which  are  considered  in  the  present  study,  the  rate 
effect  was  indeed  observed1101  but  it  is  assumed  to  have  a  limited  influence  on 
the  extrapolation  of  the  fundamental  phenomena  as  observed  in  the  model,  to 
larger  scale  systems.  Similarly,  tAlrapolation  is  not  likely  to  be  affected  by  the 
size  effect  law  considerations  when  the  structure  failure  is  dominated  by 
yielding  of  the  reinforcing  steel.  Undoubtedly,  further  study  of  these 
phenomena  will  deepen  the  understanding  and  knowledge  of  their  influence  on 
both  the  soil-structure  problem,  and  on  the  model  to  prototype  extrapolation 
of  other  experimental  systems  of  similar  nature. 
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Therefore,  the  value  of  the  results  from  the  small  scale  experimental  system  of  the 
present  study,  lies  mainly  in  revealing  the  mechanisms  that  are  involved  in  its 
response  and  their  influence,  as  discussed  in  the  previous  sections.  The  above 
discussion  shows  that  although  specific  extrapolation  might  not  be  practical  (and  was 
indeed  not  intended),  these  mechanisms  as  analyzed  here  are  likely  to  be  involved 
in  similar  larger  scale  soil-embedded  structure  systems. 


Table  2-1 

Properties  of  the  Micro-Concrete  Slabs 


iul  radius,  R  »  2.50  inches. 

100-A,/sd,  in  each  direction,  where  A,  »  9.5115  10"4  in1,  is  the  cross  sectional  area  of  the  #20  wire  that  was  used,  a  is  the  spacing,  and  d  *  h  -  0.06  inch. 
31000  psi. 


Table  2-2 

Sequence  of  the  Impactor  Hits 
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Table  2-3 

Method  of  Data  Acquisition 

Information  /  Record 

Source 

Start  time  of  the  vertical  velocity  and 

Top  accelerometer  output  (singly  and 

displacement  at  roof  center. 

doubly  integrated,  respectively). 

Amplitude  of  the  vertical  velocity  and 

LVDT  output  (first  derivative  and 

displacement  at  the  roof  center. 

direct  output,  respectively;  time  is 

shifted  according  to  the  above). 

Start  time  of  the  whole  structural 

Bottom  accelerometer  output  (once 

vertical  velocity  and  displacement. 

and  double  integrated,  respectively). 

Note:  Since  the  magnitudes  of  the  rigid  body  velocity  and  displacement  are  less  than  10%  of 
the  roof  center  velocity  and  displacement  (compare  test  results  in  figures  2-1 5  and  2-16),  the  difference 
between  the  relative  and  absolute  motion  records  of  the  roof  center  is  not  large.  However,  for 
determining  the  possible  load-displacement  (or  load-velocity)  relationship  at  the  roof  edge,  the  rigid 
body  displacement  record  (bottom  accelerometer)  is  not  negligible,  and  even  if  its  amplitude  had 
relatively  larger  experimental  errors  (than  the  LVDFs  records),  important  information  about  the  start 
time  of  that  displacement  and  its  general  shape  (in  the  time  domain)  is  provided. 
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Table  2-4 

Wave  Velocities  and  Propagation  Time  in  Micro-Concrete  Slabs 


. 

E  =  2-106  psi 

E  =  S-lO^si 

v  =  0.1 

1 

B 

HDEH 

v  =  0.2 

cL  (in/sec) 

102419 

106752 

161938 

168790 

tr  (h  =  0.4") 

7.8  Msec 

7.5  Msec 

4.9  Msec 

4.7  Msec 

tr  (h  =  0.5”) 

9.7  Msec 

93  Msec 

6.2  Msec 

5.9  Msec 
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Table  2-5 

Extrapolation  Ratios  from  a  Small  Scale  System  to  a  Larger  Scale  System 


Description  of  the  variable 

Dimension 

Prototype/Model 

ratio 

Geometry  variables  which 
have  a  length  dimension 

L 

n 

Geometry  variables  which 
have  area  dimension 

L2 

n2 

Variables  which  have 
dimension  of  force  per 
area  (e.g.  pressure,  bulk 
modulus). 

FL'2 

1 

Non-dimensional  variables 
(e.g.  strain,  Poisson  ratio) 

— 

1 

Variables  which  have 
dimension  of  mass  density 

ML  3  = 
VT2La 

1 

Accelerations 

LT2 

n'1* 

Velocities 

LT1 

1 

Time 

T 

n 

•  with  the  exception  of  the  acceleration  of  gravity  (i.e.,  ^  =  gp). 


ftxperimenl ;»l  scl  up  for  a  small  scale  model  of  a  soil  embedded  structure  system 
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Fig.  2-6:  Setup  for  measurement  of  the  free  field 
vertical  stress  distribution. 
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I  pressure;  DOB  =  2.5  inches,  H  =  25  inches. 


74 


Percent  of  Microaggregate 


Sieve  # 

30 

50 

100 

100 


Percent  sand 
by  weight 

20 

70 

10 


Gradation  Chart  of  the  Sand 
Sieve  Number -U.S.  Standard  Series 


Fig.  2-9:  Sieved  portions  of  the  sand  used  for  the  micro-concrete  specimens. 
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Fig.  2-10-a:  Form  and  reinforcement  ready  for  casting  of  micro-concrete  slabs. 


w 


Fig.  2-10-b:  Casting  of  micro-concrete  slabs  -detailed  view  of  the  double 
reinforcement  placed  in  the  form  (The  spot  welded  meshes  are  held 
apart  with  vertical  spacers). 


Fig.  2-10-c:  Removal  of  the  specimen  from  the  form  after  casting. 


Fig.  2-10-d:  Casting  of  micro-concrete  slabs  -  the  final  product,  a  doubly  reinforced 
micro-concrete  slab  (Actual  dimensions:  0.45  inch  thick,  6  inches 
diameter). 


hydraulic 


11:  Setup  for  testing  the  micro-concrete  slabs  under  direct  uniform  pressure. 


Fig.  2- 12- a:  Details  of  a  structure  -with  a  SIMPLY  SUPPORTED  roof. 


Fig.  2-12-b:  Detail  of  the  model  structure  with  a  SIMPLY  SUPPORTED  roof. 
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Fig.  2-13-b:  Detail  of  the  model  structure  with  a  CLAMPED  roof. 


0.0025 


Free  Vibration  of  a  simply  supported  slab. 
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Fig.  2-14-b:  Free  Vibration  of  a  clamped  slab. 
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Pressure  and  velocity  signals  of  the  plexiglas  roof. 
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Fig.  2-19:  Impulse  records  over  the  plexiglas  roof. 
(DOB  =  3  inches,  H  =  25  inches) 


(0.&4%) 


center  displacement;  inch 

Load-displacement  curves  of  micro-concrete  plates  under 
pressure  applied  at  different  rates. 


milliseconds 

Pressure  induced  under  the  target  plate  center 
n  the  tests  with  DOB  =  3  inches  (H  =  25  inches). 
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Fig.  2  -29:  Simply  supported  v.s.  clamped  roof  -  impulse  records. 


milliseconds 

-30:  Impulse  records  over  roof  slabs  at  different  DOB. 


CHAPTER  3 


NUMERICAL  ANALYSIS 

3.1  Description  of  the  Numerical  Method. 

Numerical  analysis  of  the  system  was  performed  by  employing  a  nonlinear,  two 
dimensional  dynamic  finite  element  code,  SAMSON2[36J.  The  program  was  run  on 
an  HP9000  series  300  Work-station,  on  an  "Encore  Multimax  320  machine",  and  on 
a  "CRAY-2"  super-computer.  The  numerical  simulation  serves  as  an  analytical  tool, 
supplementary  to  the  experimental  system.  The  axisymmetric  geometry  of  the 
experimental  setup  allows  a  numerical  analysis  with  a  two  dimensional  model.  The 
main  features  in  the  SAMSON2  code  which  were  used  for  this  part  of  the  research 
are  the  following: 

1.  Explicit  time  integration,  using  a  central  difference  method.  The  program 
automatically  evaluates  a  maximum  time  step  size  for  a  stable  solution.  This  value 
is  usually  an  upper  bound  because  it  includes  neither  damping  nor  nonlinear  effects. 

2.  Nonlinear  material  modeling  (with  a  tension  cutoff),  which  was  used  for  the  soil 
and  for  the  reinforced  concrete  roof.  The  material  models  are  discussed  separately 
in  the  following  sections. 
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3.  Slide  lines  for  modeling  material  interfaces.  The  slideline  is  an  interface 
defined  by  the  user.  It  separates  two  groups  of  elements,  allowing  each  group  to 
have  different  element  size,  such  that  the  number  of  nodes  on  each  side  of  the 
slideline  is  not  necessarily  the  same.  There  are  two  possible  motion  relations 
between  the  nodes  of  each  side  of  the  slideline:  A  tied  slideline  forces  one  set  of 
nodes,  defined  as  "slave  nodes",  to  have  displacements  which  are  dictated  by  the  set 
of  nodes  located  on  the  other  side  of  the  slideline,  and  defined  as  "master  nodes". 
The  other  possibility  is  an  interface  such  as  in  the  soil-structure  case  along  which,  slip 
or  separation  may  occur.  Both  types  of  slidelines  were  employed  for  this  analysis  as 
can  be  seen  in  the  finite  element  meshes  (e.g.,  figures  3-7  and  3-11).  Although  the 
slide  line  interface  is  convenient,  it  does,  however,  require  input  of  a  "capture 
distance",  to  determine  if  nodes  on  opposite  sides  of  an  interface  are  close  enough 
to  act  together.  This  numerical  parameter  has  been  observed  to  influence  the 
response  of  the  system.  It  should  be  noted  that  during  the  calibration  process  of  the 
code  (described  in  the  following  sections)  a  single  value  of  capture  distance 
parameter  was  used.  Thus,  the  numerical  results  are  suitable  for  the  analysis  of  the 
mechanisms  of  the  soil-structure  system  but  should  be  considered  with  judgement 
when  and  if  used  for  design. 

Further  improvement  of  the  numerical  slideline  interface  feature,  although  beyond 
the  scope  of  this  research,  can  be  useful  for  design  and  analysis  of  soil-structure 
systems  which  indeed  include  such  interfaces. 


The  numerical  study  included  three  parts: 


Calibration  of  the  material  models:  The  experimental  results  with  the  elastic-linear 
Plexiglas  roof,  with  the  reinforced  micro-concrete  roof,  and  of  the  free  field  (see 
Chapter  2)  were  used  as  references  for  the  calibration  of  the  material  models. 
Numerical  simulations  of  these  experiments  assisted  in  the  determination  of  the 
material  parameters  for  the  soil,  the  structure,  and  for  the  interface  between  them. 

Supplementary  research:  The  data  acquisition  system  and  motion  and  load 
transducers  provided  rather  detailed  experimental  results  (Chapter  2).  Nevertheless, 
some  information,  such  as  the  shear  stresses  in  the  soil  above  the  structure,  or  the 
pressure  in  the  soil  around  the  structure  was  not,  or  could  not  be  measured.  The 
numerical  analysis  supplemented  the  experimental  research  by  providing  these  data 
as  well  as  the  necessary  information  to  confirm  some  of  the  doubts  which  arose 
concerning  the  interpretation  of  the  accelerometers  readings  of  the  absolute  motion 
of  the  roof  and  of  the  whole  structure. 

Larger  scale  systems:  The  program  was  used  to  model  with  a  plane  geometry  a 
larger  scale  experiment.  The  analysis,  described  in  Chapter  4,  was  performed  in 
order  to  examine  the  physical  mechanisms,  which  were  observed  in  the  small  scale 
system  (Chapter  2).  While  the  geometry  and  the  material  parameters  of  the  larger 
scale  experiment  were  different  from  those  which  were  used  in  the  small  scale 


system,  similar  behavior  was  noted. 
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33,  Numerical  Models  and  Simulations. 

3 3  1  External  Input  Load  for  the  Soil-Structure  System. 

The  external  load  was  introduced  through  10  "load-lines".  Each  "load-line” 
describes  the  pressure  as  a  piece-wise  linear  function  of  time.  The  pressure  under 
the  target  plate  was  measured  in  the  experiments  described  in  section  2.2.1.  With 
the  assumption  of  an  axisymmetric  spatial  distribution,  the  loaded  surface  was  divided 
into  annular  strips  according  to  the  location  of  the  surface  nodes  in  the  finite  element 
mesh  (Figure  3-1).  Over  each  annular  strip,  i.e.  between  every  two  loaded  nodes,  the 
pressure  signal  was  divided  into  several  time  steps,  thus  simulating  the  actual 
measured  load. 

3.2.2  Numerical  Model  of  Plexiglas. 

The  plexiglas  structure  was  modeled  by  axisymmetric  conical  shell  elements.  These 
elements  can  be  used  in  the  program  only  with  an  elastic-plastic  biaxial  material 
model.  According  to  the  assumption  of  plane  stress  the  stresses  normal  to  the  shell 
plane  are  zero.  The  following  stress-strain  relationships  for  the  shell  elements  relate 
to  figure  3-2: 

E 

o,_  =  - 

“  (1-v2) 


<e«+veee) 


(3.1) 
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°ee  =  — ^T<vejM+eee)  (3*2) 

(1-v2) 

Ojy  =  0  (3.3) 

trr  =  “  (T^)(e«+e6fl)  (3'4) 

The  following  material  properties  are  provided  by  the  manufacturer  for  plexiglas: 

density  p  =  1.12  104  lb-se^/in4 

Young’s  modulus  E  =  3.5  ‘10s  -  5  10s  psi 
Poisson’s  ratio  v  =  0.35 

Results  of  the  numerical  simulation  of  the  free  vibration  of  a  simply  supported 
Plexiglas  slab  were  compared  with  experimental  measurements.  The  Plexiglas  plate 
was  modeled  by  five  axisymmetric  shell  elements  loaded  by  a  short  duration  load  at 
the  center.  A  similar  frequency  of  the  mid-point  vibration  was  obtained  with  an 
elastic  modulus  of  3.25 ’105  psi  (figure  3-3). 


The  disadvantage  of  this  material  model  is  its  inability  to  model  accurately 
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materials  that  are  weak  in  tension  (such  as  soil  or  concrete).  This  disadvantage 
becomes  insignificant  when  the  model  is  used  for  the  elastic  Plexiglas  structural 
elements.  The  biaxial  state  of  stress  in  the  shell  elements  is  another  disadvantage, 
which  becomes  important  when  a  transverse  tension  wave  is  likely  to  propagate 
through  the  thickness  of  the  roof  slab  (see  also  2.4.1).  This  disadvantage  does  not 
apply  to  the  main  part  of  the  numerical  analysis,  since  continuum  elements  were  used 
to  model  the  concrete  roof  rather  than  type  of  material  model  and  element 
mentioned  above. 

3 23  Numerical  Model  of  Sand. 

Sand  is  a  granular  material,  which  when  dry,  has  little  or  no  tension  capacity.  It 
is  a  medium  which  when  displaced,  particularly  under  dynamic  conditions,  is  not 
necessarily  a  continuum.  Therefore,  the  geometrical  modeling  of  sand  by  finite 
elements  is  a  numerical  constraint  of  the  actual  granular  material.  Figure  3-4 
illustrates  the  assumptions  fundamental  to  the  numerical  modeling  of  the  soil. 

The  sand  was  modeled  by  4-node  quadrilateral  continuum  elements  with  one 
integration  point.  Its  mass  density  was  measured  and  found  to  be  equal  to  1.6  10^ 
lb-sec2/ in*.  The  material  model  for  the  sand  was  the  "AFWL  Engineering  Model" 
given  by  this  code,  and  corrected  and  modified  for  this  research.  The  user  prescribes 
hydrostatic,  non-linear,  stress-strain  relations  and  a  yield  surface,  which  may  be 
enhanced  according  to  the  strain  rate.  The  sign  convention  specifies  compression 
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stresses  and  strains  as  negative.  Under  applied  loading  (i.e.,  in  compression,  the 
volumetric  strain  is  negative,  and  its  absolute  value  increases),  the  volumetric  stress 
is  determined  according  to  the  volumetric  strain  and  the  prescribed  loading  curve. 
The  position  of  the  unloading  curve  is  updated  according  to  the  minimal  strain  and 
stress  during  its  loading  history  and  therefore  shifted  horizontally,  parallel  to  the 
volumetric  strain  axis  and  retains  the  tensile  cutoff  at  its  prescribed  level.  For  re¬ 
loading,  the  volumetric  stress-strain  follows  the  unloading  curve.  When  and  if  the 
volumetric  strain  exceeds  the  minimum  level  ever  reached  during  loading,  the 
volumetric  stress  is  determined  again  from  the  loading  curve.  The  parameters  that 
need  to  be  quantified  are: 

Constitutive  relations: 


K  = 


for  loading, 
for  unloading. 


(3.5) 


for  loading. 
v(<g  for  unloading. 


(3.6) 


Note  that  K(cv)  for  unloading,  may  have  a  "cutoff  value";  i.e.  K(crv>avcu,off)  =  0 
means  no  strength  in  tension,  beyond  a  volumetric  stress  higher  than  avcu,off. 
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The  failure  surface  (which  may  be  enhanced  according  to  the  strain  rate)  is: 

y  =  =  n- oj  (3.7) 

In  equations  3.5,  3.6,  and  3.7: 

K  -  bulk  modulus 

€v  -  volumetric  strain  (positive  for  expansion) 
av  -  normal  pressure  (positive  for  tension) 
v  -  Poisson’s  ratio 
Y  -  failure  surface 

(J2’)1/2  -  second  deviatoric  stress  invariant 

For  an  axisymmetric  geometry,  a  cylindrical  coordinate  system  is  used,  and  the 
stresses  increments  for  each  time  step  are: 

Aou  =  2p(Aeu  -  -^)  +  Aov  ;  i  =  r,  B,  y  (3  8) 


A°r,  * 

2  (1-v) 


(3.9) 


where: 


is  the  shear  modulus. 
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K  is  the  volumetric  stress  increment. 


Ac  =  Ac  +Ae.„+Ac  is  the  volumetric  strain  increment. 

v  rr  80  yy 

At  each  time  step  these  incremental  stresses  are  calculated  according  to  the 
volumetric  strain  and  stress,  and  the  current  material  parameters  (equations  3.5, 3.6, 
3.8  and  3.9).  The  second  deviatoric  stress  invariant  is  determined  by: 

V^=/^(orr-oee)2+(oe0-°,/+(oyy'orr)2l+ory  (3,10) 

and  is  compared  to  the  corresponding  failure  surface  level,  Y(-av).  If  exceeded,  then 
the  stresses  are  modified  to  reduce  the  second  deviatoric  stress  invariant  to  the 
failure  surface  level  while  maintaining  the  same  volumetric  stress: 

Y 

oti  »  ov  +  —  (ou  -  oj  ;  i  =  r,  0,  y  (3.11) 


°n  = 


- 

& 


(3.12) 
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323.1  Literature  survey  of  experiments  with  20-30  Ottawa  Sand. 

For  a  first  approximation  some  experimental  data  with  20-30  Ottawa  sand  were 
reviewed.  Then  the  free  field  experimental  results  (described  in  section  2.2.2)  were 
used  for  the  determination  of  the  sand  modeling.  The  following  experimental  results 
were  gathered  and  reviewed  in  order  to  establish  some  typical  constitutive 
relationships  as  well  as  the  failure  envelope  for  Ottawa  sand: 

Ko  and  Scott1291,  tested  dense,  medium-dense,  and  loose  sand  under  cyclic 
hydrostatic  compression.  Their  results  provide  the  constitutive  relationship  and  give 
the  change  of  volumetric  strain  as  a  function  of  volumetric  stress.  Since  the  tests 
were  performed  with  cyclic  loading,  they  provide  the  necessary  information  for 
determining  the  bulk  modulus,  K,  during  loading  and  unloading  (as  a  function  of  the 
volumetric  strain  and  the  volumetric  stress,  respectively),  for  loose,  medium-dense, 
and  dense  Ottawa  sand. 

More  experimental  data  were  provided  from  tests  done  by  Ueng,  Tzou,  and  Lee1501. 
They  made  triaxial  tests  with  Ottawa  sand.  Their  reported  results  were  digitized  and 
transformed  into  plots  of  volumetric  strain  and  volumetric  stress.  Schmertmann1381 
made  CFS  (Cohesion-Friction-Strain)  tests,  which  indicate  for  dense  saturated 
Ottawa  Sand  under  axial  compressive  strain  larger  than  0.5%,  cohesion  and  a  friction 
angle  of  about  0.05  Kg/cm2  (0.71  psi)  and  33  degrees,  respectively.  They  also 
provided  the  failure  condition  for  the  confining  pressure  at  which  they  were 
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conducted  (2  Kg/cm2,  28.44  psi).  More  data  for  the  failure  conditions  (i.e.  failure 
envelope)  were  taken  from  the  experimental  results  reported  by  Ko  and  Scott1281,  for 
dense,  medium-dense,  and  loose  Ottawa  sand.  The  experiments  were  performed 
under  static  conditions.  Selig  and  Vey1391  tested  Ottawa  sand  under  static  and 
dynamic  induced  stresses.  While  the  purpose  of  the  dynamic  tests  was  to  study  stress 
propagation  and  attenuation,  the  static  tests  provided  the  failure  conditions  for  three 
confining  pressures,  5,  7.5,  and  12.5  psi.  It  was  reported  that  in  some  cases  the 
dynamic  peak  pressure  was  raised  above  the  static  strength,  while  for  other  cases  it 
was  lower.  Whitman  and  Healy1531  also  examined  the  response  of  sand  to  dynamic 
loading  and  reported  "little,  if  any,  strain-rate  effect  on  the  strength  of  dry  Ottawa 
Sand". 

3232  Calibration  of  the  numerical  model  for  sand. 

The  free  field  stress  distribution  3  inches  under  the  target  plate  was  used  as  a 
reference  for  the  material  parameters  of  the  model.  Figures  3-5  through  3-8 
summarize  the  experimental  results  with  Ottawa  sand.  They  also  include  the  actual 
parameters  that  were  used  to  calibrate  the  numerical  model  according  to  the  free 
field  experimental  results.  The  finite  element  mesh  which  was  used  (shown  in  figure 
3-9)  describes  half  of  the  experimental  set-up  (applying  axial  symmetry).  It  included 
216  nodes,  176  elements,  and  2  slidelines.  The  upper  part  of  the  mesh  included 
elements  which  modeled  "loose"  sand,  while  the  other  elements  modeled  "dense" 
sand.  The  external  load  simulated  the  experimental  measurements  of  the  impact 
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under  the  target  plate  as  described  in  section  3.2.1.  A  constant  Poisson’s  ratio  of  0.2;> 
was  assumed. 

Figures  3-10  through  3-12  describe  the  soil  pressure  ob;ained  numerically  versus 
the  experimental  results,  under  an  impact  generated  by  a  drop  of  25  inches.  The 
vertical  pressure  is  compared  at  3  inches  under  the  target  plate  at  locations  at  the 
center  line,  2  inches  from  the  center,  and  2.5  inches  from  the  center.  The  soil 
parameters  yielded  numerical  predictions  in  good  agreement  with  the  experimental 
measurements  and  they  were  used  in  the  next  steps  of  the  FEM  analysis. 

3.2.4  Numerical  Simulation  of  the  Soil-Elastic  Structure  System. 

The  next  step  examined  the  model  containing  a  buried  structure.  Figure  3-13 
describes  the  mesh  which  was  used  to  simulate  the  experiment  which  was  performed 
with  the  plexiglas  roof  slab  (section  2.2.3).  It  had  282  nodes,  181  elements,  and  8 
slideline  interfaces.  The  Plexiglas  structure  was  modeled  with  axisymmetric  shell 
elements,  section  3.2.2..  Vertical  stresses  over  the  roof  as  well  as  the  motion 
(displacement  and  velocity)  of  the  roof  and  of  the  whole  structure,  were  computed 
by  the  finite  element  program,  modeling  the  small  scale  experiments  described 
previously  (2.2.3).  Among  other  parameters  of  the  system  examined  was  the  soil-roof 
slideline  interface.  The  following  comparison  was  performed  between  the  numerical 
and  the  experimental  records:  Displacement  of  the  center  of  the  roof,  pressure  at 
the  center  of  the  roof,  and  pressure  over  its  edge,  2.5  inches  from  the  center.  The 
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input  was  the  load  function  used  for  the  free  field  runs.  Modeling  of  the  sand  was 
taken  according  to  results  obtained  from  the  free  field  runs,  figures  3-14  through  3- 
17. 

The  numerical  calculations  model  well  the  experimental  results  up  to  about  1  msec. 
Then  the  roof  center  rebounded  sharply  upwards  in  the  numerical  model,  while  it 
maintained  a  lower  frequency  response  in  the  experiment  (figure  3-14).  This 
behavior  was  reflected  in  the  velocity  record  and  also  caused  a  sharp  increase  in  the 
pressure  records  over  the  structure  (see  FEM  records  at  t>l  msec  in  figures  3-15,  3- 
16,  and  3-17,  respectively). 

The  runs  with  the  Plexiglas  roof  combined  the  soil  and  the  structure  and  provided 
optimal  calibration  of  the  capture  distance  value  of  the  slide  line  interface  between 
them.  No  friction  or  adhesion  was  introduced  at  this  interface,  and  the  value  of  the 
capture  distance  was  kept  unchanged  in  the  following  runs  with  the  reinforced 
concrete  slab  as  the  structure  roof.  These  runs  also  provided  additional  information 
to  that  which  was  obtained  experimentally.  As  it  can  be  seen  in  figure  3-14  the 
absolute  center  displacement  and  the  displacement  relative  to  the  floor  are  similar. 
Hence,  the  records  of  the  relative  motion,  experimentally  obtained  by  the  LVDT,  are 
verified  *s  representative  to  the  motion  of  the  roof  center  since  the  relatively  smaller 
absolute  motion  of  the  whole  structure  was  indicated  by  the  accelerometers  signals 
(see  Chapter  2)  as  well  as  in  the  numerical  results.  Further  numerical  results  of  the 
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soil-Plexiglas  roof  run  are  discussed  in  section  3.3. 

32.5  Numerical  Model  of  Reinforced  Concrete. 

No  attempt  was  made  to  use  a  sophisticated  material  model  for  the  concrete.  The 
material  model  which  was  used  merely  applied  the  bas'c  behavior  of  concrete,  and 
the  analysis  was  focused  on  the  soil-structure  interaction  mechanisms  which  are 
involved.  Modeling  of  the  reinforced  micro-concrete  slabs  was  performed  with  the 
4-node  quadrilateral  continuum  elements  for  concrete,  and  with  shell  elements  for 
steel  (figure  3-18).  While  concrete  is  continuous  and  can  be  geometrically  modeled 
by  continuum  elements,  it  is  necessary  to  use  an  equivalent  geometrical  model  for 
the  reinforcing  steel  meshes.  The  axisymmetric  shell  elements  represent  a  continuous 
disk  equivalent  to  the  discrete  steel  mesh,  where  the  equivalence  was  based  on  equal 
force  taken  by  the  numerical  and  the  actual  steel.  Considering  a  cross  section  of  the 
slab  at  an  arbitrary  angle  o,  the  internal  total  force  taken  by  the  steel  bars  over  the 
cross  section  is  Fx  and  Fy  in  the  x  and  y  directions,  respectively  (Figure  3-19).  The 
number  of  steel  bars  in  each  direction  over  the  cross  section,  Nx  and  Ny  depends  on 
the  spacing,  s.  For  an  equal  spacing  in  each  direction,  they  are: 

^  2lsina 

*  *  (3.13) 

..  2/cos a 
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For  reinforcement  bars  of  cross  sectional  area  A,,  and  assuming  equal  stress  in  the 
x  and  y  direction  (4x=fsy=fs),  the  total  forces  Fx  and  Fy  are: 

F  -N  A J  =  ~/2lsma 

s  (3.14) 

Fy=NyA/s  =  ^/,2/COSa 


The  total  force  normal  to  the  cross  section,  Fn,  should  be  equal  to  the  force  taken 
by  the  equivalent  steel  disk  with  a  thickness  h^: 

Fn=Fx  sina  +  Fy  cosa  =  2//r^  (3.15) 


Substituting  Fx  and  Fy  from  equation  3.14  into  equation  3.15,  and  arranging  terms, 
the  equivalent  thickness  for  the  reinforcement  shell  elements  is  obtained: 


(3.16) 


The  material  model  for  the  reinforcing  shell  elements  was  the  biaxial  elasto-plastic 
model  with  an  assumed  elastic  modulus  of  30 -106  psi,  and  the  measured  yield  stress, 
fy,  of  51000  psi  (see  also  chapter  2).  It  should  be  noted  that  the  shell  elements  were 
too  thin  to  contribute  significant  flexural  stiffness. 
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Because  of  the  non-linear  nature  of  the  constitutive  relations  of  concrete,  and  its 
relatively  low  strength  in  tension,  it  was  also  modeled  by  the  "AFWL  Engineering 
Model"  (described  in  section  3.2.3)  with  suitable  parameters  for  concrete. 
Calibration  of  this  material  model  was  performed  with  initial  reference  to  the 
hydrostatic  relations  and  failure  surface  based  on  Chen191  and  Ottosen[35L  The  yield 
surface  normalized  with  respect  to  the  compressive  strength,  fc  is  based  on  the 
following  function: 
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(3.17) 


where:  (J2)1/2  is  the  second  deviatoric  stress  invariant. 

Ij  is  the  first  stress  invariant;  Ij  =  3 av 
av  is  the  volumetric  stress. 
fc  is  the  uniaxial  compressive  strength, 
a,  /3,  y,  &  are  parameters  which  are  fitted  with  experimental  results. 

The  parameter  y  =  0  at  failure  and  the  other  parameters  which  are  suggested  for 
failure  conditions  in  the  literature  are  1.3,  9,  and  3.2  for  a,  0,  and  6,  respectively. 

The  curves  which  were  used  for  the  analysis  are  shown  in  figures  3-20  and  3-21. 
They  were  evaluated  by  two  separate  comparisons  with  experimental  records  of 
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micro-concrete  slabs:  The  free  vibration  of  the  simply  supported  slab  of  test  7  (see 
Table  2-1)  was  recorded  (before  the  test  was  done).  The  slab  was  numerically 
modeled  by  a  finite  element  mesh  as  described  in  figure  3-18  with  the  dimensions  of 
the  slab  in  test  No.  7.  The  comparison  between  the  numerical  and  experimental 
record  is  shown  in  figure  3-22.  Because  of  the  small  displacement  (about  0.002 
inches)  the  range  of  the  strains  and  stresses  in  the  slab  was  relatively  low  (maximum 
compressive  strain  of  less  than  0.001),  and  far  from  failure.  Thus,  it  provided  the 
initial  parameters  for  the  material  model.  Poisson’s  ratio  was  assumed  to  be  equal 
to  0.2  and  constant.  The  initial  elastic  modulus  of  micro-concrete  was  —  1 -106  psi, 
section  2.2.4.2.  Figure  3-22  also  shows  the  effect  of  the  lateral  restraint  of  the 
support  on  the  response  of  the  slab  (dashed  line  in  figure  3-22).  As  it  can  be  seen 
the  natural  frequency  is  doubled  when  the  support  is  hinged  but  its  lateral 
displacement  is  restrained. 

These  model  parameters  were  then  completed  and  refined  by  using  the  high  rate 
experimental  results  of  the  experiments  with  circular  slabs  that  were  described  in 
section  2.  The  loading  pressure  and  the  experimental  versus  numerical  record  of  the 
mid-point  displacement  under  this  pressure  are  given  in  figures  3-23  and  3-24, 
respectively.  In  figure  3-24  the  numerical  modeling  of  the  reinforced  micro-concrete 
slab  with  the  above  parameters  for  the  "AFWL  Engineering  Model",  was  seen  to  be 
in  good  agreement  with  the  experimental  results.  It  should  be  noted,  however,  that 
once  the  concrete  continuum  is  cracked,  its  finite  element  representation  becomes 
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less  accurate.  Although  the  material  law  provides  a  tension  cut-off,  a  tension  failure 
which  forms  a  physical  crack  cannot  be  accurately  modelled.  Therefore,  an 
equivalent  average  value  was  used  for  the  tension  cut-off  that  was  lower  than  the 
actual  tension  capacity  of  the  concrete.  Damage  or  failure  condition  is  evaluated  by 
the  strain  at  the  reinforcing  steel  elements  and  by  comparison  of  the  slab 
displacement  to  a  relevant  load-displacement  curve,  but  the  extent  of  the  damage  (or 
failure)  at  very  large  displacements  (  >  1/10  of  the  span),  is  not  modeled.  This 
solution  was  found  to  be  satisfactory  for  the  soil-structure  modeling  provided  the 
response  of  the  embedded  structure  depends  mainly  on  early  time  mechanisms,  which 
may  cause  later  larger  displacements.  For  the  current  soil-structure  system  it  is, 
indeed,  the  early  time  response  which  is  more  important,  and  therefore  the  above 
model  is  adequate  for  its  analysis. 

3 2.6  Numerical  Simulation  of  the  Soil-Micro-concrete  Structure  System. 

The  element  types  which  are  described  in  the  previous  sections  were  used  for  the 
numerical  simulation  of  test  No.  6  (see  Table  2-1),  with  corresponding  material 
models.  The  input  pressure  applied  at  the  surface  under  the  target  plate  was 
generated  under  a  height  of  drop  of  25  inches.  The  experimental  and  numerical 
input  under  the  center  of  the  target  plate  are  shown  in  figure  3-25.  Figure  3-26 
shows  the  finite  element  mesh  which  was  employed  for  this  analysis.  The  analytical 
results  are  plotted  with  their  corresponding  experimental  records  in  figures  3-27 
through  3-30.  In  the  early-time  response  there  is  a  good  agreement  between  the 
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analytical  and  experimental  results.  The  center  displacement  and  velocity  follow 
their  corresponding  experimental  records  during  the  whole  significant  part  of  the 
response  (figures  3-27, 3-28).  As  in  the  case  of  the  elastic  Plexiglas  roof,  the  relative 
and  absolute  displacements  are  similar,  indicating  that  the  displacement  of  the  whole 
structure  was  much  smaller  compared  to  the  displacement  of  the  roof  center.  The 
velocity  record  (figure  3-28)  also  indicates  possible  instability  of  the  numerical 
solution  at  t  >  0.6  msec,  as  the  signal  becomes  relatively  noisy.  Such  a  problem  might 
be  addressed  by  further  decrease  of  the  time  steps.  However,  it  was  not  done 
because  the  time  step  size  which  was  used  was  already  small  ( 1 -1CT7  sec)  and 
required  a  long  computation  time.  Additionally,  the  essential  features  of  the  early- 
time  part  of  the  response  was  already  captured  by  the  numerical  solution. 
Accordingly,  the  pressure  in  the  soil  element  above  the  roof  agrees  well  with  the 
corresponding  experimental  records  up  to  t«0.6  msec.  Although  high  peaks  are 
shown  which  were  not  measured  in  the  experiment,  they  became  zero  at 
approximately  the  same  time  as  in  the  experiment  (-1.2  msec  in  figures  3-29,  and  3- 
30). 
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33  Discussion  of  the  Numerical  Simulations  for  the  Experimental  Results. 

The  FEM  analysis  of  the  small  scale  system  experiments  enabled  the  parameters 
of  the  code  to  be  calibrated,  and  also  provided  additional  information,  difficult  to 
acquire  experimentally.  Records  of  vertical  pressure  and  shear  stresses  in  the  soil, 
together  with  the  vertical,  transverse  stresses  through  the  micro-concrete  slab  mid¬ 
section,  enable  further  examination  of  the  mechanisms  that  influence  its  response 
(see  Chapter  2). 

33.1  Displacement  of  the  Roof  Slab. 

The  numerical  displacement  solution,  which  follows  the  experimental  measurement 
during  about  2  msec,  also  indicates  that  the  early-time  response  is  significant,  and 
that  the  mechanisms  which  take  place  during  that  time,  up  to  -0.6  msec  in  the 
current  system,  are  dominant  in  influencing  the  total  response.  These  mechanisms 
are  described  in  Chapter  2.  At  t=0.4  msec  (figure  3-14)  the  central  displacement 
was  ~0.02  in  and  the  strain  at  the  reinforcing  steel  element  at  the  bottom  center  of 
the  slab  exceeded  the  yielding  strain,  0.0017  (E  is  assumed  30-106  psi  for  steel,  and 
fy  =  51000  psi).  The  radial  strain  and  stress  in  the  bottom  steel  element  were  0.0021 
and  51000  psi,  respectively.  At  t=0.6  msec  the  displacement  at  the  center  was  0.036 
in  and  the  radial  strain  in  the  bottom  steel  element  was  0.0053.  Referring  to  the 
load-displacement  curves  of  similar  slabs  (figure  2-23),  the  yielding  of  the  bottom 
reinforcing  steel  element  at  a  central  displacement  of  0.02  in  is  reasonable,  and 
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indicates  initiation  of  flexural  damage. 

3.3.2  Transverse  Stresses  through  the  Slab. 

Figure  3-31  shows  the  transverse  stresses  in  the  slab  mid-section,  top,  middle,  and 
bottom  elements,  during  the  first  0.4  msec  of  the  slab  response.  It  can  be  seen  that 
as  the  incident  pressure  wave  propagated  through  the  thickness  of  the  slab  it 
attenuated  and  reached  a  minimum  level  at  the  bottom  element.  Considering  the 
fact  that  the  stresses  in  the  quadrilateral  elements  are  calculated  at  their  one  mid¬ 
integration  point,  it  is  understood  that  there  is  some  non-zero  vertical  stress  in  the 
slab  bottom  element,  and  that  its  relatively  low  value  agrees  well  with  the  natural 
boundary  condition  of  zero  vertical  stresses  at  the  free  bottom  of  the  slab.  This 
boundary  condition  can  be  satisfied  by  a  reflection  tension  wave  from  the  free 
bottom  of  the  slab. 

As  the  slab  center  starts  to  gain  velocity,  just  before  t  =  0.2  msec,  the  vertical  stress 
at  the  top  of  the  slab  reaches  its  peak  and  begins  to  drop,  while  at  the  same  time 
maintaining  a  steady  minimum  level  at  the  bottom  (figure  3-31).  The  stresses  at  the 
soil  above  the  slab  also  start  to  decrease  at  the  same  time  (figure  3-29). 

333  Soil-Structure  Displacements. 

One  of  the  purposes  of  using  the  slideline  interface  feature  in  the  numerical  model 
was  to  examine  the  relative  motion  of  the  soil-structure.  However,  due  to  the 


120 


incorporation  of  the  capture  distance  parameter,  this  part  of  the  analysis,  might  not 
be  reliable.  Figure  3-32  describes  the  displacement  of  two  nodes  0.5  inches  from  the 
center,  located  on  opposite  sides  of  the  sand-roof  interface.  There  appear  to  be 
some  gaps  between  the  "sand  node"  and  the  "roof  node";  however  they  are  of  the 
order  of  10-4  inch.  This  difference  is  more  likely  to  be  a  numerical  gap  rather  than 
to  represent  a  physical  one.  Therefore,  as  in  the  experimental  observations,  a  zero 
stress  was  assumed  at  the  soil  element  above  the  structure  to  represent  a  possible 
physical  gap. 

3.3.4  Shear  Stresses  in  the  Soil  above  the  Structure. 

Figures  3-33  and  3-35  show  the  early  time  shear  stresses  in  the  soil  at  the  level 
above  the  Plexiglas  and  the  micro-concrete  roofs,  respectively.  Note  that  the  clear 
span  of  the  Plexiglas  roof  was  2.75  inches  while  that  of  the  micro-concrete  roof  was 
2.5  inches.  From  these  figures  it  can  be  seen  that  at  a  time  of  0.3  to  0.4  msec 
relatively  high  shear  stresses  developed  on  the  sides  of  the  structure.  Comparing 
figures  3-33  and  3-35  with  figures  3-16  and  3-29,  that  time  is  seen  to  be  when  the 
pressure  over  the  center  of  the  roof  was  decreasing  and  reached  its  lower  zero  value. 
In  the  case  of  the  micro-concrete  roof  the  pressure  over  the  roof  edge  peaked  at  0.3- 
0.4  msec  (figure  3-30)  and  over  the  Plexiglas  roof  it  peaked  for  the  second  time  after 
previously  having  decreased  (figure  3-17).  These  observations  support  the 
interpretation  of  the  experimental  results,  that  arching  develops  in  the  soil  over  the 
structure  after  the  first  peak  of  the  pressure  has  already  loaded  the  structure.  The 
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shear  stresses  over  the  wall  of  the  structure  coincide  with  an  arch  action  at  which 
support  there  is  a  relatively  high  shear.  (In  the  asymmetrical  case  there  is  actually 
a  dome  action  analogous  to  the  arch  action  which  in  the  following  text  is  referred  to 
as  arch  action.).  In  the  soil  above  the  walls  of  the  structure  there  is  a  transition  from 
a  relatively  low  level  of  vertical  stresses  to  a  zone  of  high  vertical  stresses,  either 
above,  or  adjacent  to  the  stiffer  support  of  the  structure.  The  location  of  these  zones 
might  differ  according  to  the  stiffness  of  the  roof  and  to  the  foundation  stiffness, 
which  affects  the  deflection  of  the  whole  structure.  The  change  of  sign  of  the  shear 
stresses  in  this  zone  is  associated  with  a  lateral  displacement  of  the  sand  located 
within  the  high  pressure  zone,  on  both  sides  of  the  roof.  From  figures  3-33  and  3-35 
it  can  be  seen  that  while  this  region  was  located  on  both  sides  of  the  Plexiglas 
structure,  -3  inches  from  the  center  (in  figure  3-33),  in  the  case  of  the  stiffer 
concrete  roof  it  was  above  the  walls,  -2.5  inches  from  the  center  (in  figure  3-35). 

The  arch,  or  dome  action  in  the  soil  at  t  >  0.35  msec  is  also  demonstrated  in  figures 
3-34  and  3-36,  which  show  the  vertical  pressure  distribution  above  the  Plexiglas  and 
the  concrete  roof  level,  respectively.  A  concentration  of  vertical  stresses  in  the  soil 
above  the  relatively  stiffer  wall  area  indicates  the  formation  of  arch  action  over  the 
deflecting  center  area  of  the  structure. 

3 3.5  Effect  of  the  Boundary  Conditions. 

The  case  of  test  No.  6  in  which  the  roof  was  simply  supported,  was  simulated  with 
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a  clamped  type  of  support:  Additional  Plexiglas  shell  elements  attached  to  the 
concrete  element  at  the  bottom  edge  and  to  the  top  Plexiglas  wall  element  increased 
the  bending  stiffness  at  the  support.  Figure  3-37  shows  the  center  displacement  of 
the  roof  in  the  two  cases.  The  clamped  roof  (solid  lines  in  figure  3-37)  reached  a 
peak  deflection  of  about  85%  of  that  of  the  simply  supported  roof  at  about  the  same 
time.  Although  the  smaller  displacement  of  the  clamped  roof  is  expected,  the 
relatively  small  difference  in  the  deflection  indicates  either  that  the  clamping  model 
was  not  very  effective  or  that  plastic  hinges  developed  at  the  slab  edges  at  an  early 
time  of  the  response.  In  either  case  the  soil  pressure  records  over  the  roof  (figures 
3-38  and  3-39)  show  the  same  response  during  the  first  peak  with  higher  pressure 
over  the  clamped  roof  center  in  the  later  time  part  of  the  response.  This  observation 
indicates  an  influence  of  the  boundary  conditions  on  the  arching  action  in  the  soil 
over  the  structure.  However,  because  arching  takes  place  at  a  later-time  of  the 
response  (after  the  first  peak),  its  effect  is  significant  only  if  the  structure  sustains  the 
early-time  peak  loading  by  causing  further  reduction  of  the  load. 

3 3.6  Effect  of  the  Soil-Structure  Interface  Properties. 

In  the  numerical  analysis,  which  was  described  in  the  previous  sections,  sliding  was 
allowed  along  the  soil-structure  interface.  A  similar  run  was  performed  with  fixity 
at  the  soil-structure  interface.  The  results  are  shown  in  figure  3-40  through  3-42. 
The  effect  of  the  properties  of  the  interface  is  displayed  by  a  smaller  deflection  of 
the  roof  and  by  corresponding  smaller  load  over  the  roof  center  in  the  fixed  case 


(solid  line  versus  dashed  line  in  figures  3-40  and  3-41).  It  can  be  seen  in  figure  3-41 
that  the  load  over  the  center  of  the  roof  increases  to  the  same  peak  in  both  cases, 
but  drops  faster  in  the  case  of  a  fixed  interface.  At  t>0.3  msec  the  load  over  the 
edge  in  the  fixed  case  is  higher  than  that  in  the  case  of  free  sliding  interface  (solid 
versus  dashed  lines  in  figure  3-42).  These  observations  indicate  an  earlier 
development  of  soil  arching  when  the  soil-structure  interface  is  less  compliant.  The 
arch  action  is  caused  by  the  deflection  of  the  structure,  and  by  the  soil  shear  capacity 
and  is  enhanced  by  increasing  the  friction  between  the  soil  and  the  structure.  It 
should  be  noted,  though,  that  this  indication  must  be  further  investigated  for  two 
main  reasons: 

1.  While  the  soil-structure  interface  friction  may  allow  earlier  development  of 
arching  in  the  soil,  it  also  loads  the  structure  with  shear  stresses  at  its  interface.  The 
above  numerical  result  should  be  followed  by  further  experimental  study  of  the 
structure  response  to  increased  interface  shear  stresses. 

2.  Since  the  interface  model  involves  a  numerical  parameter  of  capture  distance 
this  result  indicates  a  qualitative  trend,  which  must  be  verified  and  quantified  by 
further  experimental  research  and  improvement  of  the  numerical  modeling  of  the 


soil-structure  interface. 


Numerical  modeling  of  soil. 
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note:  unloading  curve  is  adjusted  by 
horizontal  shift  according 
to  the  volumetric  stress. 
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Volumetric  Strain  (compression)  , 

Loading  and  unloading  curves  for  numerical  modeling  of  Ottawa  sand. 


Volumetric  Pressure  (positive  for  compression),  psi 

Fig.  3-8:  Yield  surface  for  Ottawa  sand  - 
experimental  measurements  and  actual  values  used  for  the  FEM  analysis 
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Fig.  3-10:  Free  field  pressure  at  the  center  -  > 
numerical  v.s.  experimental  results  (DOB  =  3  inches;  H  =  25  inches). 
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Fig.  3-14:  Displacement  of  the  roof  center  -  FEM  analysis  v.s. 
experimental  results  (Plexiglas  roof;  DOB  =  3  inches;  H  =  25  inches) 
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15:  Velocity  of  the  roof  center  -  FEM  analysis  v.s.  experimental  results 
(Plexiglas  roof;  DOB  =  3  inches,  H  =  25  inches) 
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18:  Finite  elements  mesh  model  of  a  simply  supported  circular  reinforced  concrete  slab. 
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(Volumetric  Pressure  /  f'c)  ;  (positive  for  compression) 

Fig.  3-21:  Yield  surface  for  the  numerical  modeling  of  micro-concrete. 
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Vibration  of  a  simply  supported  micro-concrete  slab. 
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3-24:  Mid-point  displacement  of  a  micro-concrete  slab  under  a  high  rate 
uniformly  distributed  pressure  (  f’c  =  6500  psi,  v  -  0.2,  h  =  0.47  inch) 
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-26:  Finite  element  mesh  for  modeling  of  test  No. 
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FEM  analysis  v.s.  experimental  results  (test  6,  2nd  drop): 
Displacement  of  the  roof  center 
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FEM  analysis  v.s.  experimental  results 
Velocity  of  the  roof  center. 
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-29:  FEM  analysis  v.s.  experimental  results  (test  6,  2nd  drop). 
Pressure  at  the  CENTER  of  the  roof 
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analysis  v.s.  experimental  results  (test  6,  2  drop) 
Pressure  over  the  EDGE  of  the  roof 
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Fig.  3-32:  Vertical  displacements  at  the  sand-roof  interface, 
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Fie.  3-34:  Pressure  distribution  above  the  rool 
umerical  analysis;  Plexgiglas  roof;  DOB  =  3  inches,  H 
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.  3-36:  Pressure  distribution  above  the  roof  level, 
umerical  analysis  of  m.c.  roof,  test  6,  2“d  hit) 
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Numerical  analysis  of  a  simply  supported  v.s.  clamped  micro-concrete 
roof  -  vertical  pressure  at  the  CENTER  of  the  roof. 
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-39:  Numerical  analysis  of  a  simply  supported  v.s.  clamped  micro-concrete 
roof  -  vertical  pressure  over  the  EDGE  of  the  roof. 
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CHAPTER  4 


SIMULATION  AND  ANALYSIS  OF  A  LARGER  SCALE  SYSTEM 


The  finite  element  program,  which  is  described  in  Chapter  3,  was  employed  for  the 
numerical  analysis  of  a  soil-embedded  structure  system  of  different  geometry,  larger 
scale  and  with  higher  external  load. 

4.1  The  Reference  Svstem  and  ..s  Numerical  Model. 

4.1.1  Description  of  the  Soil-Structure  System. 

For  this  purpose,  an  experiment  conducted  by  Kiger,  Eagles  and  Baylot1271  was  used 
as  a  reference.  They  performed  a  series  of  tests  with  a  massive  post-tensioned 
concrete  reaction  structure  that  clamped  its  roof  slab  structural  elements  for  rotation 
as  well  as  lateral  displacement  and  was  designed  to  ensure  one  way  action  of  the  test 
slabs.  The  reinforced  concrete  slabs  were  2.9  inches  thick,  and  had  a  clear  span  of 
24  inches  together  with  a  0.5%  tension  and  compression  reinforcement  with  0.5  inch 
cover.  Test  number  4  of  this  series  was  performed  with  a  dynamic  external  load  and 
sand  backfill,  and  is  numerically  simulated  here.  The  structure  was  buried  at  a 
shallow  DOB  equal  to  half  of  the  roof  span  (i.e.,  12  inches),  and  subjected  to  an 
external  high  explosive  blast  simulating  a  nuclear  weapon.  The  experiment  report 
has  the  following  information  for  the  materials: 
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1.  Concrete  was  reported  to  have  an  average  28-days  compressive  strength  of  6900 
psi  and  modulus  of  elasticity  equal  to  4.56 -lO6  psi. 

2.  Reinforcement  was  made  of  small-diameter  deformed  wire,  0.177  inch  diameter 
with  an  average  yield  strength  of  90200  psi. 

3.  Backfill  of  test  number  4  was  a  locally  available  sand  called  Reid-Bedford 
model  sand.  Its  average  mass  density  was  1.53  -10-4  lb-sec2/in4. 

4.1.2  Numerical  Modeling  of  the  Soil-Structure  System. 

4.1.2.1  Geometrical  modeling. 

Figure  4-1  describes  the  geometrical  modeling  of  the  reference  system  by  a  finite 
element  mesh,  which  represents  half  of  the  system  by  symmetry  and  assuming  plane 
strain.  To  simplify  the  calculations,  only  the  roof  slab  is  modeled,  while  the 
supporting  structure  is  represented  by  lateral  and  vertical  restraint  at  the  structure 
bottom  edge  node  and  a  lateral  restraint  at  its  upper  edge  node.  Plane  continuum 
elements  were  used  for  modeling  the  sand  and  the  concrete;  and  bar  elements  were 
used  to  model  the  reinforcement.  The  area  of  the  bar  element  was  determined  by 
analysis  similar  to  that  described  in  Chapter  3.  Requiring  equal  force  to  be  taken 
by  the  numerical  reinforcement  model  and  by  the  actual  reinforcing  mesh,  leads  to 
an  equivalent  numerical  reinforcement  area  per  bar  element,  A^: 
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A  =  =  Ai*a*ai'f*  =  p-1  ”>d  =  — •  1-2.44  =  0.012—  (4.1) 

**  fs  ft  100  * 


4.122  Material  modeling. 

In  this  example  the  materials  which  are  included  in  the  system,  reinforced  concrete 
and  sand,  are  similar  to  those  of  the  smaller  scale  system.  Therefore,  the  material 
law  which  was  used  to  model  the  sand  and  the  concrete  was  the  "AFWL  Engineering 
Model",  described  in  the  previous  chapter.  The  steel  was  modeled  by  a  uniaxial 
stress  material  law  with  parameters  suitable  to  the  actual  steel  wires  used  in  the 
experiment  (see  section  4.1.1).  The  parameters  used  for  the  "AFWL  Engineering 
Model"  are  shown  in  figures  4-2  to  4-5  and  are  not  the  same  parameters  which  were 
used  for  the  small  scale  system  but  were  evaluated  as  follows: 

1.  The  load  on  the  system  was  much  higher  than  the  low  impact  generated  by  the 
drop  of  the  impactor  steel  ball  (peak  of  3300  psi  versus  600  psi,  respectively). 
Therefore,  the  constitutive  relationships  in  the  high  stress  range  of  the  soil  had  to  be 
evaluated  on  the  basis  of  typical  stress-strain  relationship  of  granular  soils,  as 
suggested  by  Lambe  and  Whitman1311  (see  figures  4-4  and  4-5):  For  stresses  up  to 
about  2000  psi,  the  stress-strain  curves  are  concave  upward,  as  the  soil  undergoes 
what  is  called  "locking".  Starting  at  about  2000  psi,  the  stress-strain  curve  begins  to 
develop  a  reverse  curvature  and  becomes  concave  to  the  strain  axis.  During  this 
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stage,  sand  particles  fracture  and  permit  tighter  packing  which  again  makes  the 
stress-strain  curve  concave  at  high  stresses  and  strains.  Fracturing  of  the  sand 
particles  becomes  important  when  the  stresses  exceed  5000  psi. 

The  initial  parameters  of  the  curve  are  similar  to  those  shown  in  figure  3-5  for 
loose  Ottawa  sand  and  were  based  on  the  arrival  time  of  the  pressure  wave  as 
reported  in  the  actual  test.  Poisson’s  ratio  was  assumed  to  be  0.25  (as  for  the  Ottawa 
sand). 

2.  The  concrete  parameters  shown  in  figures  4-2  and  4-3  were  based  on  the  same 
principles  at  those  used  for  the  modeling  of  the  micro-concrete19351  with  the  reported 
initial  Young’s  modulus  of  4.56-106  psi.  The  failure  surface  (figure  4-3)  was  based 
on  the  reported  compressive  strength  of  6900  psi.  Poisson’s  ratio  was  assumed  to  be 
0.2. 

4.1 23  Load  modeling. 

Because  of  the  nature  of  the  explosive  charge  used,  the  load  was  assumed  to  be 
uniformly  distributed.  The  pressure  signal,  as  recorded  and  reported,  was  digitized 
(from  Ref.  27)  and  used  as  an  input  at  the  top  surface  of  the  finite  element  mesh. 
Figure  4-6  shows  the  time  record  of  the  input  pressure. 
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4*2  Results* 

The  results  of  the  numerical  analysis  which  was  performed  with  the  parameters 
described  in  section  4.1.2  are  discussed  here.  They  are  compared  with  the 
experimental  records  which  were  reported  in  Ref.  27  and  further  analyzed  in  view 
of  the  conclusions  drawn  from  the  small  scale  system  experiments  (Chapters  2  and 
3). 

4.2.1  Comparison  of  the  Numerical  and  Experimental  Results. 

Figure  4-7  shows  the  displacement  record  of  the  mid-roof.  Both  the  numerical  and 
experimental  records  indicate  failure  of  the  roof  slab.  The  numerical  curve  is 
different  from  the  experimental  in  the  time  of  its  initial  response  and  in  its  maximum 
deflection.  While  the  numerical  analysis  shows  a  mid-point  deflection  which  starts 
at  t«0.5  msec,  the  experimental  curve  shows  no  apparent  deflection  up  to  about  1 
msec  later  (i.e.  t— 1.5  msec).  Both  reach  a  1  inch  deflection  within  about  2.5  msec 
(from  the  time  deflection  started),  but  while  the  experimental  gage  measured  a 
maximum  deflection  of  about  3  inches,  the  numerical  record  showed  a  maximum 
displacement  of  about  1  inch.  These  differences  should  be  judged  in  view  of  the 
notes  which  were  reported  (Ref.  27):  The  gage  which  was  used  for  the  displacement 
measurements  was  reported  to  be  insensitive  to  deflections  occurring  in  the 
millisecond  range  which  explains  the  difference  in  the  response  time  between  the 
numerical  and  the  experimental  record  (see  also  in  Appendix  B  a  similar  observation 
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about  the  LVDT  gage  in  the  small  scale  system).  Furthermore,  it  was  reported  that 
in  the  experiment  the  gage  probe  bent  and  the  gage  failed.  Moreover,  failure  at  very 
large  displacements  is  not  modelled  well  in  the  analysis,  as  noted  in  section  3.2.5. 
Direct  loading  of  a  similar  slab  was  tested  prior  to  the  dynamic  experiment1271  and 
reported  to  cause  failure  at  a  center  displacement  of  about  0.25  in.  The  load- 
deflection  curve  which  was  reported  (insert  in  figure  4-7)  was  obtained  under  static 
loading.  This  indicates  that  although  the  numerical  model  may  not  predict  the  actual 
total  displacement  at  failure,  it  shows  a  deflection  higher  than  0.25  in,  accompanied 
by  yielding  of  the  steel  elements,  which  implies  flexural  damage  (cMeci>cyiel(j  at 
t=0.78  msec  and  at  t=2.5  msec  at  the  support  and  the  center,  respectively). 

The  experimental  normal  pressure  at  the  soil-structure  interface  is  compared  to  the 
vertical  stress  in  the  elements  above  the  roof  in  figures  4-8  to  4-11.  The  pressure 
gage  which  was  located  at  the  roof  center  showed  an  initial  spike  in  front  of  the  main 
peak  of  the  pressure  signal  (not  shown  in  figure  4-8)  and  the  gage  was  reportedly 
broken  very  early  in  the  test.  Thus  the  earlier  numerical  arrival  time  at  the  center 
is  reasonable,  and  its  peak  agrees  well  with  the  experimental  one.  Over  the  rest  of 
the  slab  roof  there  is  reasonably  good  agreement  between  the  numerical  and  the 
experimental  pressure  records  during  the  main  first  peak  of  the  loading  (figures  4-9, 
4-10,  and  4-11).  The  relatively  higher  secondary  peak  in  the  numerical  result  may 
be  due  to  the  difference  between  the  actual  supporting  structure  and  its  numerical 
modeling.  While  the  actual  structure  was  embedded  in  the  backfill  and  probably 
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underwent  some  rigid  body  displacement  (not  recorded  in  the  experiment),  it  was 
modelled  in  the  finite  element  mesh  by  a  rigid  support  (figure  4-1).  Hence,  possible 
later  time  displacement  of  the  whole  structure  might  have  caused  arching  in  the  sand 
above  it  which  reduced  the  total  loading  over  the  supports  of  the  deflecting  structure. 
In  the  model,  however,  such  displacement  was  not  possible,  therefore  the  higher 
vertical  pressure  over  the  rigid  support  during  1.5  msec<t<3  msec. 

422  Discussion. 

As  in  the  case  of  the  small  scale  system,  the  experimental  results  are  limited  to  the 
number  and  capability  of  the  measuring  equipment.  The  numerical  finite  element 
analysis,  after  careful  calibration,  can  provide  further  information  usually  difficult  to 
obtain  in  the  experiment. 

Figure  4-12  shows  the  vertical  transverse  stresses  across  the  mid-section  of  the  roof 
together  with  the  soil  vertical  pressure  above  it.  While  at  the  interface  the  vertical 
stress  is  similar  to  the  soil  pressure,  it  drops  to  a  minimum  level  at  the  bottom  of  the 
slab  where  the  natural  boundary  condition  requires  zero  normal  pressure.  This  is 
similar  to  the  numerical  results  of  the  small  scale  system  (see  also  3.3.1).  At  t-0.6 
msec  the  roof  center  starts  to  gain  velocity  (figure  4-13)  and  the  pressure  in  the  soil 
above  it  starts  to  drop.  This  pressure-velocity  early  time  dependency  was  observed 
also  in  the  small  scale  system  (Chapters  2  and  3)  and  seems  to  dictate,  together  with 
the  free  field  pressure  level,  the  subsequent  total  response  of  the  roof. 
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The  numerical  stress  distribution  over  the  roof  is  plotted  in  figure  4-14.  It  is 
interesting  to  see  that  the  initial  pressure  peaks  uniformly  over  the  entire  roof  while 
it  is  lower  in  the  soil  adjacent  to  the  roof  (t=0.7  msec  in  figure  4-14  and  t=0.65  msec 
in  figure  4-15).  Then,  within  only  0.1  msec  the  normal  pressure  drops  sharply  over 
the  roof  center,  while  the  drop  over  the  support  (located  12  inches  from  the  center) 
is  much  milder,  and  in  fact  remains  higher  than  the  free  field  level  in  the  soil 
adjacent  to  the  structure  (at  a  distance  of  15  inches  from  the  center  in  figure  4-14 
and  at  t~l  msec  in  figure  4-15).  The  initial  increase  of  the  pressure  over  the  roof  is 
typical  of  a  rigid  boundary  effect  on  an  incident  pressure  wave.  Although  the 
pressure  dropped  rapidly,  this  initial  increase  of  the  pressure  was  not  observed  in  the 
small  scale  system,  and  the  difference  between  the  initial  response  of  the  two  systems 
can  be  explained  by  the  difference  in  the  external  load.  The  high  explosive  charge 
in  the  larger  scale  experiment  generated  a  higher  level  of  pressure  than  the  small 
scale  system  and  also  had  a  shorter  rise  time  to  its  peak.  Therefore,  while  the  initial 
load  level  in  the  small  scale  system  was  still  relatively  low  when  the  structure  began 
to  respond,  it  was  significantly  higher  in  the  larger  scale  system. 

In  spite  of  this  difference,  the  rest  of  the  response  resembles  that  of  the  small  scale 
system.  This  is  demonstrated  in  figure  4-15  where  the  numerical  results  for  the 
pressure  and  motion  of  the  roof  were  plotted  with  the  experimentally  reported  free 
field  pressure  (assumed  uniformly  distributed).  Points  "A"  and  "B",  which  were 
defined  by  the  observations  of  the  small  scale  system  results  (Chapter  2,  section  2.4), 
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are  also  identified  here:  As  the  velocity  of  the  roof  center  develops,  the  center  load 
decreases  and  keeps  dropping  while  the  pressure  over  the  support  increases  at  t=0.8 
msec  (point  "B")  above  the  reference  free  field  level,  indicating  the  development  of 
arching.  More  evidence  for  this  mechanism  is  shown  in  figure  4-15,  where  the  shear 
stresses  in  the  soil  above  the  structure  are  plotted.  Similar  to  the  results  of  the  small 
scale  system,  the  arching  action  is  indicated  here  by  the  development  of  relatively 
high  shear  stresses  in  the  soil  element  over  the  supports  (see  also  section  3.3.3). 

The  development  of  arching  in  this  case,  at  t*l  msec,  occurred  when  cracking  had 
already  developed  near  the  support  and  started  to  develop  at  the  center  (center 
deflection  of  -0.2  in  in  figure  4-7),  and  is  another  example  of  the  conclusions  that 
can  be  drawn  from  the  small  scale  experiment  (Chapter  2)  that  the  early-time 
response  in  a  soil-embedded  structure  systems  is  dominant.  The  zero  to  peak  portion 
of  the  free  field  incident  wave  at  the  structure  level  hits  the  structure,  and  is  even 
amplified  until  a  tension  wave  is  reflected  from  the  bottom  free  surface  of  the 
structural  element  (within  about  0.01  to  0.7  msec  in  concrete  slabs  up  to  100  in 
thick).  If  the  structure  can  sustain  this  portion  of  the  load,  then  the  rest  of  the 
mechanisms  which  are  involved  in  the  system  are  important  for  further  release  of  the 
load  off  the  structure,  or  maintaining  a  mild  degree  of  damage  in  a  worse  case. 


LOADED  SURFACE 


Fig.  4-1:  Finite  element  mesh  for  the  numerical  modeling  of  the  test 
which  was  performed  by  Kiger,  Eagles,  and  Baylot!*7' 
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Fig.  4-3:  Numerical  modeling  of  a  larger  scale  test 
Yield  surface  of  concrete. 


note:  unloading  curve  is  adjusted  by 
horizontal  shift  according 
to  the  volumetric  stress. 


0.015  0.020  0.02L  0.0J0  0.035  0.040  0 

Volumetric  Straii  (compression) 

Numerics!  modeling  of  a  larger  scale  test(r#f  ,00) 
Loading  and  unloading  curves  of  the  sand. 


Volumetric  Pressure  (positive  for  compression 

4-5:  Numerical  modeling  of  a  larger  scale  test<r,f  l00)  - 
Yield  surface  of  the  sand. 


3400 


180 


Tsd 

« 


o 

O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

O 

o 

o 

O 

O 

o 

O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

CM 

o 

00 

CD 

CM 

O 

CO 

CD 

CM 

o 

oo 

CD 

M- 

CM 

K) 

K> 

CM 

CM 

CM 

CM 

CM 

V— 

milliseconds  ,  ,  k 

Fig.  4-6:  Numerical  modeling  of  a  larger  scale  test/  11  - 
Input  pressure  at  the  surface  (assumed  uniformly  distributed). 
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Fig.  4-12:  Numerical  modeling  of  a  larger  scale  test127’  - 
Vertical  (transverse)  stresses  across  the  roof  center  cross-section. 
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Fig.  4-16:  Numerical  modeling  of  a  larger  scale  test2  - 
Shear  stresses  in  the  soil  above  the  roof. 


CHAPTER  5 


CONCLUSIONS 


5*1  Response  of  A  Soil-Structure  System, 

The  experimental  and  numerical  study  of  the  response  of  a  structure  embedded  in 
a  shallow  depth  of  burial  to  an  external  surface  impact  was  described  in  previous 
chapters.  From  the  investigation  of  the  physical  mechanisms  which  influence  the 
behavior  of  the  system,  two  main  phases  of  the  response  were  observed:  an  early- 
time  phase  and  a  later-time  phase.  The  early-time  response  of  the  system  is  usually 
during  the  first  peak  of  the  incident  impact,  which  is  first  amplified  at  the  initially 
rigid  soil-structure  interface,  and  later  attenuated  with  the  arrival  of  a  tension  wave 
reflected  from  the  free  bottom  surface  of  the  roof.  Further  decrease  of  the  initial 
load  is  due  to  the  interaction  between  the  loading  soil  and  the  loaded  structure, 
which  is  affected  first  by  the  velocity  of  the  structure  and  then  by  its  deflection.  The 
later-time  phase  of  the  response  involves  possible  gap  between  the  soil  and  the 
center  portion  of  the  roof,  which  releases  the  load  over  that  region,  and  rigid  body 
settlement  of  the  structure  which  tends  to  decrease  the  arching  action  at  shallow 
DOB,  increase  the  load  over  the  roof  center,  and  decrease  it  over  the  supports. 

According  to  the  peak  pressure  level  and  the  rise  time  of  the  incident  impact, 
whether  the  structure  will  sustain  the  load  or  fail  is  determined  already  during  the 
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first  early-time  phase  of  the  response.  The  initial  performance  of  the  structure 
depends  on  its  capacity  to  sustain  impact  loading  and  on  both  the  level  of  the  load 
and  on  its  rise  time  and  duration.  Therefore,  for  very  short  rise  times  (about  1/4  to 
1/2  of  the  natural  period)  the  total  impulse  might  need  to  be  considered  as  a  loading 
criterion  for  failure  (rather  than  just  the  peak  pressure). 

The  propagation  speed  of  transient  elastic  waves  in  the  concrete  roof  is  relatively 
high,  because  of  the  large  stiffness  of  concrete.  Since  the  thickness  of  the  roof  is, 
usually,  relatively  small,  the  reflection  of  a  tension  wave  from  the  bottom  free  surface 
of  the  roof  occurs  soon  after  it  is  hit  by  the  pressure  wave  (order  of  magnitude  of 
0.010  msec).  When  it  reaches  the  upper  surface  of  the  roof  the  amplified  incident 
pressure  is  released  approximately  to  the  level  of  the  free  field  pressure.  The  tension 
in  the  lower  surface  of  the  roof  might  also  cause  spalling  of  concrete,  which  should 
be  considered  during  the  design  phase. 

The  subsequent  interaction  mechanisms  are  important  in  decreasing  the  load  and 
thereby  limiting  the  impulse  acting  on  the  structure.  The  general  trend  of  this 
portion  of  the  response  is  to  release  the  load  from  the  centrally  loaded  zone  of  the 
structure,  and  to  lower  the  pressure  over  it  relative  to  the  free  field  pressure.  The 
drop  of  the  pressure,  which  is  faster  than  the  free  field  drop,  is  initiated  by  the 
downward  movement  of  the  structure  and  corresponds  to  its  velocity.  Further  drop 
of  loading  over  the  central  zone  of  the  deflecting  roof  is  caused  by  arching,  which 
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develops  in  non-cohesive  soils  over  the  relatively  flexible  displacing  structure.  The 
release  of  the  pressure  over  the  center  of  the  roof  due  to  arching  is  typically 
accompanied  by  a  corresponding  increase  of  the  load  over  its  stiffer  supports.  In  the 
absence  of  these  mechanisms  the  structure  would  have  been  loaded  over  a  longer 
time,  and  probably  experienced  further  damage. 

It  was  also  observed  that  the  soil-structure  interface  friction  influences  the 
response.  More  interface  friction  enhances  the  arching  mechanism  and  causes 
further  decrease  of  the  load  (after  thA  first  peak)  over  the  center  region  of  the  roof. 
However,  interface  friction  increases  the  shear  loading  at  the  structure  upper  surface, 
which  should  also  be  considered. 

The  boundary  conditions  of  the  roof  may  contribute  to  its  capacity  in  sustaining  the 
initial  portion  of  the  loading.  They  are  less  important  in  the  later-time  response, 
especially  when  multiple  hits  are  considered.  In  the  case  of  a  repeated  loading, 
plastic  hinges  at  the  clamping  supports  produce  a  response  which  is  similar  to  the 
response  of  a  non-clamping  support,  as  the  deflection  shape  of  both  become  similar. 

Finally,  the  rigid  body  settlement  motion  of  the  whole  structure  needs 
consideration,  as  it  has  a  similar  effect  on  the  later-time  velocity-load,  and  arching 
mechanisms.  Although  the  motion  develops  later  than  that  of  the  roof,  it  can  further 
contribute  to  the  process  of  load  release  over  the  supports. 
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Application  of  the  Small  Scale  System  Observations  and  of  the  Computer  Code. 

This  work  employed  an  experimental  portion  which  was  performed  in  the 
laboratory  with  a  small  scale  model  system.  The  associated  dimensional  analysis 
shows  the  ratios  for  an  extrapolation  of  the  small  scale  results  into  a  larger  scale 
system.  It  was  shown  that  if  the  gravitational  portion  of  the  involved  mechanisms  is 
not  dominant,  the  distortion  in  such  an  extrapolation  is  limited  to  the  order  of 
magnitude  of  the  experimental  error. 

The  axisymmetric  system  which  was  used  did  not  model  a  specific  larger  scale 
prototype.  However,  when  compared  with  another  larger  scale  system  of  the  same 
nature,  similar  mechanisms  were  observed.  The  extrapolation  ratios  together  with 
the  expression  for  the  amount  of  distortion  (due  to  gravitational  effects)  propose  a 
useful  way  of  using  the  results  of  a  small  scale  model  system  for  a  larger, 
geometrically  similar  system. 

One  of  the  aspects  of  the  small  scale  modeling  was  the  use  of  micro-concrete 
reinforced  with  steel  wires  to  model  reinforced  concrete.  The  employment  of  micro¬ 
concrete  in  this  work  was  based  on  previous  comprehensive  research  which 
investigated  the  optimal  way  to  model  reinforced  concrete  in  small  scale  systems.  It 
was  further  based  on  static  tests  with  micro-concrete  slabs  (which  were  performed 
during  this  research).  The  micro-concrete  showed  general  mechanical  behavior 
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which  is  typical  of  concrete.  It  should  be  noted  that  there  were  indications  that  the 
elastic  modulus  of  the  micro-concrete  was  lower  than  that  of  concrete.  Therefore, 
it  is  recommended  to  examine  and  verify  adequately  the  initial  elastic  modulus  of  the 
modeling  material  together  with  the  testing  of  its  compressive  strength. 

The  "SAMSON2"  computer  code  which  was  developed  in  the  Air-Force,  offers 
flexible  features  for  geometrical  and  material  modeling.  The  main  features  which 
correspond  to  modeling  of  soil-structure  systems,  are  the  material  models  and  the 
sliding  interfaces.  However,  these  features  are  not  developed  yet  to  their  full 
potential.  The  sliding  interfaces,  which  are  important  in  the  modeling  of  the  contact 
surface  between  two  different  materials,  should  be  developed  further  and  improved 
to  enable  better  definition  and  determination  of  the  "capture  distance"  parameter, 
to  limit  its  influence  on  the  results. 


APPENDIX  A 


MEASUREMENTS  OF  LOADS 


Two  commercially  available  load  transducers  were  considered  for  the 
measurements  of  the  pressure  induced  by  the  sand  on  the  roof  of  the  structure: 

1.  "Kulite"  pressure  transducer,  was  successfully  used  in  the  Foam  HEST 
experiments11819,251,  and  was  designed  to  measure  pressure  at  the  soil-structure 
interface.  The  major  drawback  of  this  transducer  with  regard  to  the  current 
experimental  setup  is  its  1  inch  diameter  size  and  that  it  must  also  be  mounted 
within  the  roof  of  the  structure  so  that  its  face  is  leveled  with  the  roof  upper 
surface.  For  these  reasons  it  was  used  only  as  a  reference  transducer  for 
another,  smaller,  load  cell  which  was  modified  adequately  to  match  the 
manufacturer’s  instructions. 

2.  "Entran"  ELF-500- 100  load  cell,  was  successfully  used  by  Chen[81  to  measure 
the  load  over  the  center  of  the  roof  in  a  former  stage  of  this  experimental 
program.  This  load  cell  is  only  0.5  inch  in  diameter  and  0.11  inch  in  thickness, 
has  low  mass  and  it  is  capable  of  measuring  dynamic  load  signals  in  excess  of  20 
KHz.  According  to  the  manufacturer’s  instructions1171,  for  use  over  parts  of  the 
roof  which  are  off  center,  all  nonuniform  and  horizontal  components  of  the 
pressure  must  be  isolated  to  insure  that  the  load  cell  is  only  loaded  by  a  uniform 
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normal  pressure  over  its  load  button. 

In  a  series  of  tests  several  techniques  were  examined  in  order  to  obtain  a  reliable 
pressure  signal.  The  criterion  for  the  quality  of  the  load  cell  output  was  a  reference 
signal  from  a  "Kulite"  pressure  gauge  (Figure  A-l),  which  was  acquired  by  using  a 
small  "thumb  tack"  type  of  steel  support  glued  to  the  loading  button  of  the  load  cell. 
By  mounting  the  load  cell  on  a  relatively  rigid  pin  pointed  support,  with  the  back  of 
the  load  cell  facing  the  loading  medium  (i.e.  the  sand),  the  load  which  was 
transferred  to  the  loading  button  of  the  load  cell  was  distributed  uniformly  and 
directed  normal  to  the  loading  face  of  the  transducer,  as  required. 

The  output  records  of  the  two  load  cells,  when  using  this  technique,  is  described 
in  figure  A-2.  As  it  can  be  seen  the  output  records  of  the  two  load  cells  were  very 
similar  and  was  adopted  as  the  method  with  which  pressure  was  recorded  over  the 
roof.  The  total  load  indicated  by  the  load  cell  was  divided  by  its  area  (0.1963  in2) 
to  obtain  the  average  pressure.  Additional  verification  for  the  pressure 
measurements  is  discussed  in  Chapter  2,  section  2.2. 1.1. 


APPENDIX  B 


MEASUREMENTS  OF  DYNAMIC  DISPLACEMENTS 

Two  alternative  methods  to  measure  a  dynamic  displacement  (or  its  time 
derivatives)  were  examined.  One  was  the  use  of  an  accelerometer  to  obtain  a  time 
record  of  the  acceleration,  which  can  then  be  integrated  to  give  the  velocity  and  the 
displacement.  An  alternative  way  was  the  use  of  an  LVDT  (Linear  Variable 
Displacement  Transformer)  which  gives  a  direct  record  of  the  displacement.  Time 
derivatives  of  that  record  gives  the  velocity  and  acceleration  record. 

In  order  to  examine  the  two  methods,  the  following  "calibration  test"  was 
performed:  At  the  free  end  of  an  aluminum  cantilever  beam  a  Type  4374  Bruel  & 
Kjaer  accelerometer,  and  an  LVDT  were  mounted.  The  signals  of  the  transducers 
were  recorded  by  a  4094  Nicolet  digital  oscilloscope.  In  some  of  the  tests  the  signal 
of  the  accelerometer  was  recorded  directly,  and  in  others  through  a  charge  amplifier 
(Kistler  S/N  216  model  504).  The  LVDTs  signal  was  recorded  using  a  Chaevitz 
amplifier  with  excitation  frequency  of  20  kHz  to  provide  accurate  measurements  of 
the  dynamic  signal.  The  other  end  of  the  beam  was  clamped  between  two  aluminum 
blocks  (Figure  B-l).  The  cross  section  of  the  beam  was  2  in  by  0.25  in,  and  the  span 
,  1,  was  variable. 

Neglecting  damping,  the  equation  of  motion  of  the  beam  is: 
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&  =  -±&y  (B-i) 

a*4  b 2  ar2 


where:  b2  =  -Q- 

pA 

E  -  Young’s  modulus  (llxlO6  for  aluminum). 

I  -  moment  of  inertia  of  the  cross  section. 

A  -  area  of  the  cross  section. 

p  -  mass  density  (0.093/386.4=2.407x10'*  lb-secZ/in4  for  aluminum). 

The  free  end  was  deflected  by  a  string  until  w=0.05  in  was  read  at  that  point  in 
a  dial  gage  attached  to  the  beam.  The  dial  gage  was  then  removed  and  the  string 
cut,  thus  providing  the  following  initial  conditions  : 


y(x=l,t=0)  =  w(t=0)  =  0.05  in 


(B-2) 


y(x=l,t=0)  =  w(t=0)  =  0  in/sec 


The  general  solution  for  equation  B-l  is: 
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y(x,t)  =  X(x)  (A  cos(tot)  +  B  sin(ot)  )  (B-3) 


for  which  X(x)  must  satisfy: 


m  -  k*X  =  0  (B-4) 

cbc4 


2 

where:  k4  =  — 

b2 


A  convenient  form  for  the  solution  of  equation  B-4  is: 

X(x)  =  C!(cos(kx)  +  cosh(kx))  +  C2(cos(kx)  -  cosh(kx))  + 

C3(sin(kx)  +  sinh(kx))  +  C4(sin(kx)  -  sinh(kx))  (B-5) 


The  constant  k  is  obtained  from  the  boundary  conditions.  For  a  fully  clamped 
cantilever  the  boundary  conditions  are  : 

X(0)  =  0  ;  i*(0)  =  0 


V(U)  =0  or 


=  0 

ax3 


M(lj)  =0  or 


dx 2 


From  these  boundary  conditions  we  get : 


cos(kl)«osh(kl)  =  -1  (B-6) 

The  natural  frequencies  for  the  vibration  modes  of  the  beam  are  o—blq2  (b  is 
defined  in  the  equation  of  motion,  B-l).  The  first  five  eigenvalues  for  kl  are  listed 
in  Table  B-l: 


Table  B-l 

First  Five  Eigenvalues  of  kl  for  a  Cantilever  Beam 


kjl 

k2l 

k3l 

k4l 

k5l 

1.875 

4.694 

7.855 

10.996 

14.137 

It  should  be  noted,  however,  that  a  fully  fixed  end  is  very  hard  to  obtain,  and  even 
a  slight  slope  at  the  "fixed"  support  alters  the  boundary  condition  and  actually  gives 
a  moment  of  finite  magnitude,  K  (Figure  B-l(b)).  In  that  case  the  boundary 
conditions  at  the  support  (x=0)  become: 

X(0)  =  0 

M(0,t)  =  K—( 0)  =  £7—  (0) 

a*  ax2 


While  the  conditions  at  the  free  end  remain  unchanged. 
From  these  boundary  conditions  we  get  the  following  equation: 
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(coskl +coshjfc/)^cos#  +cosh£/ +^^(-sinjW +sinhJU)  j  + 
(siniU+siiihiU)^siniU-sinhjU+-^(cosiy-coshJk/)j  =  0 


(B-7) 


The  moment  stiffness  of  the  support,  K,  can  be  either  estimated  or  found  in  a 
semi-empirical  way  by  measuring  the  frequency,  (*>,  from  which  the  value  of  k  can  be 
derived  and  substituted  into  equation  B-7.  Then,  repeating  the  test  with  a  different 
span,  1,  the  same  K  is  used  to  calculate  the  new  k,  which  corresponds  to  the  new 
frequency,  o>  for  that  span. 

Once  <i)  is  obtained,  the  solution  for  the  motion  of  the  free  end  (w  and  its 
derivatives)  is: 

w(t)  =  w0cos(<i>t) 

w(t)  =  -wwosin(a)t)  (B-8) 

w(t)  =  -<i>2w0cos(wt) 

where  w0  is  the  initial  displacement  of  the  free  end  of  the  cantilever  (0.05  inch 
in  this  case). 
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The  experiments  were  performed  with  three  spans  (I):  6.625  in,  8.0  in,  and  9.4375 
in.  From  the  output  of  the  LVDT  in  the  test  performed  with  a  span  1  =  6.625  in  (Fig. 
B-2)  the  measured  period  of  the  first  mode  was  T=  10  msec.  The  lowest  frequency 
was  then,  <t>=2ir/T= 628.32  rad/sec,  which  gives  kl  =  1.337.  Substituting  kl  into 
equation  B-7,  K= 6162.8  lb-in.  It  can  be  shown  that  this  value  of  K  allows,  for  an 
initial  deflection  of  0.05  in,  a  maximum  slope  at  the  support  of  0.8  degrees.  Such  a 
small  slope  is  indeed  undetectable  by  the  naked  eye,  but  causes  a  difference  when 
the  frequency  response  of  the  beam  is  considered.  A  fully  clamped  beam  of  the 
same  span  would  have  a  period  twice  as  long,  i.e.  Tdamped=5  msec  (see  Table  B-l). 

Using  this  value  of  K,  the  following  predictions  can  be  made  for  the  other  spans: 

For  1  =  6.625  inch  (compare  Fig.  B-2) 

Vmax  =  0.05  to  =31.4  in/sec 
amax  =  0.05  oj  2  =  51  g 

For  1  =  8.0  inch  (compare  Fig.  B-3) 
k  =  0.1729 

id  =  461.45  rad/sec,  T  =  13.62  msec. 
vmax  =  0.05  (*>  =  23.1  in/sec 
amax  =  0.05o2  =  27.55  g 
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For  1  =  9.4375  inch  (compare  Fig.  B-4) 
k  =  0.1508 

<■>  =  351  rad/sec,  T  =  17.9  msec. 

Vmax=  0.05  (o  =  17.55  in/sec 
JW  =  0.05 <i> 2  =  15.95  g 

Figures  B-2,  B-3,  and  B-4  show  the  experimental  results.  From  these  figures  the 
following  conclusions  can  be  drawn: 

1.  The  output  of  the  LVDT  is  accurate  and  stable,  and  measures  well  the 
amplitude  of  the  motion  as  well  as  its  frequency  and  can  also  be  used  to 
obtain  the  velocity  record  by  taking  the  first  time  derivative. 

2.  The  output  of  the  accelerometer,  when  amplified,  measures  well  the  frequency 
of  the  motion  signal  but  seems  to  "drift"  giving  an  unstable  integrated  record 
for  the  velocity  and  the  displacement,  after  about  half  a  period  of  the  motion. 

3.  When  not  amplified,  the  raw  signal  coming  out  of  the  accelerometer  might  be 
too  weak  (compared  to  the  noise  recorded  with  it),  thus  yielding  erroneous 
output  of  both  the  amplitude  and  the  frequency  (compare  figures  B-2,  B-3, 
with  figure  B-4). 


207 


4.  A  close  observation  of  the  two  records  shows  that  the  accelerometer  responds 
faster  than  the  LVDT  to  the  exciting  motion  by  50-100  Msec.  In  some  types 
of  experiments  (such  as  the  soil-structure  system  considered  in  this  work)  the 
relationship  between  the  motion  of  the  structure  and  the  load  acting  on  it  is 
sought,  and  this  response  time  might  be  important.  Therefore,  the 
accelerometer  should  be  used  together  with  the  LVDT  to  detect  the  actual 
starting  time  of  the  motion,  so  that  the  signal  of  the  LVDT  could  be  shifted 
accordingly. 


APPENDIX  C 


DIMENSIONAL  ANALYSIS  OF  A  SMALL  SCALE  SOIL-STRUCTURE  SYSTEM. 


Figure  C-l  describes  the  system  and  some  of  the  parameters  which  are  scaled  and 
analyzed.  Notations  of  the  dimensions  are  as  follows: 

[mass]  =  M 
[length]  =  L 
[time]  =  T 
[force]  =  F  =  M  L  /  T2 

The  variables  (and  their  dimensions)  are: 

Variables  with  length  dimension,  L: 

DOB,  h,  1,  D,  r-  see  Fig.  C-l 
u  -  average  grain  size  of  the  soil 
w  -  mid  roofs  deflection 

Variables  with  area  dimension,  L2: 

A,  -  reinforcement’s  area 

Variables  with  pressure  dimension,  F/L2  =  M/  L  T2: 


stress  in  the  soil 


or  -  reinforcement  (of  structure’s  roof)  stress 

ac  -  concrete  (of  structure’s  roof)  stress 

£  -  compressive  strength  of  concrete 

fy  -  yield  stress  of  reinforcement 

q0  -  peak  external  pressure  (see  Fig.  C-l) 

a  -  stress  acting  on  the  roof 

Ec  -  concrete  Young’s  modulus  (initial) 

Er  -  steel  Young’s  modulus 

K,  -  soil’s  bulk  modulus  for  loading 
K*  -  soil’s  bulk  modulus  for  unloading 

Variables  with  mass  density  dimension,  M/L3: 
pt  -  soil’s  density 
pc  -  roofs  (concrete)  density 

Variables  with  velocity  dimension,  L/T: 
v  -  mid  roofs  velocity  ( =  w) 

Variables  with  acceleration  dimension,  L/T2: 
g  -  gravitational  acceleration 


a  -  mid  roofs  acceleration 


Variables  with  time  dimension,  T: 


t  -  time 

to  -  typical  time  of  external  loading  (see  Fig.  C-l) 

Nondimensional  variables: 
ec  -  concrete  strain 

er  -  steel  strain 

et  -  soil  strain 

vc  -  concrete  Poisson  ratio 

vr  -  steel  Poisson  ratio 

vs  -  soil  Poisson  ratio 

f(r)  -  distribution  of  external  pressure  (see  Fig.  C-l) 


Assumptions: 


1.  All  variables  of  length  dimension  are  scaled  identically: 


=  n 


(C-l) 


where  the  subscript  "p"  denotes  prototype,  the  subscript  "m"  denotes  model,  and 
n  is  the  scale  factor  for  the  length.  Hence,  only  one  typical  length  variable,  A,  may 
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be  considered  for  the  dimensional  analysis. 

2.  Materials  of  the  prototype  and  of  the  model  systems  have  the  same  properties, 
hence: 

a.  The  stress-strain  relations  of  the  prototype  soil  is  identical  to  that  of  the  model, 
yielding: 

(^c)p  —  (®c)m»  (^r)p  —  (^r)m’  (^»)p  “(^s)m 

Note:  This  was  verified  for  the  micro-concrete  and  for  the  micro¬ 
reinforcement  in  Reference  12. 


b.  The  mass  density  of  the  soil,  concrete,  and  reinforcement  is  the  same  for  the 
model  and  for  the  prototype: 

(Pc)p  “  (Pc)m»  (Pr)p  “  (Pr)m>  (P*)p  —  (Ps)m 


Following  the  "formal  procedure"  of  a  dimensional  analysis,  the  behavior  of  the 
system  is  according  to  the  general  equation  C-2,  in  which  C  is  a  nondimensional 
constant: 


dal/-1a2_a3_a4_a5  _*6  r>a7  r  «8y  a*>r^  ctlO  ir  all  c  «12  _  «13  _  i 

<Io  °c  °t  °  *c  S  ^  fcr  P*  Pc 

yalS  ^a!6  g«17  j«18  j^«19^«20  _  £ 


a  14 


(C-2) 


Note  that  formally,  the  nondimensional  variables  may  also  be  part  of  equation  C- 
2,  but  since  they  can  not  contribute  to  the  dimensional  analysis  they  are  assumed  to 
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be  included  in  the  constant  C. 

To  keep  equation  C-2  dimensionally  consistent  the  sum  of  the  powers  of  each  of 
the  three  basic  dimensions,  L,  M,  and  T,  must  be  equal  to  zero.  Substituting  the 
dimensions  of  the  variables  into  equation  C-2  and  collecting  the  powers  of  each 
dimension,  the  three  following  equations  are  obtained: 

for  L:  atra2*Ot3*at4"O[5"Q!6‘Ot7"aE8'Q!9*O!10*fl[irQ!12*^a13*^al4‘^  *16  a17^ 

2a*  =  0  (C-3) 

for  M:  a2+a3+a4  +  as+a6+a7+a8+ci9+al0+au  +  a12+a13+a14  =  0  (C-4) 

for  T .  -2a2-2a3-2a4-2a5-2a6-2a7-2a8-2a9-2a jq-2ce j j-2qe32~® j^-2ot |^-2oc 

OE18‘^a19  =  0  (C-5) 

Any  (and  only)  three  of  the  power  a’s  may  therefore  be  determined.  Writing  the 
algebraic  relations  as  functions  of  a,,a14)  and  a18  we  get,  from  equations  C-3,  C-4, 
and  C-5: 

Qj  =  -2fl2*2ci3-2ct4-2ci5-2ci^-2ci9-2Q[8-2fl9-2otjQ-2Q!jj-2otj2-0[j5-Otig-Q[j7- 

2®20 


(C-6) 
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a14  “  ■Qt2*a3"0E4'a5‘a6"a7”0t8‘a9*at10‘Q!ir0t12‘0[13  (C-7) 


Ojg  =  2<i2 '*■20:3+20:4+ 2oi5+2ct6+2a7+ 2ctg+2oi9+2<2jo+ 2otjj  +  2ct12+ oti5  +  2o:16+ 

2ai7"ai9  (C-8) 


Take  log  of  equation  C-2  and  substitute  alt  a14,  and  a18,  to  get: 
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£  “/log71*  =  logC 

i-2 

<•14,18 


(C-9) 


The  jt  numbers  are  nondimensional,  and  provide  the  scaling  ratios  (for  a  true, 
undistorted  model)  for  each  variable,  by  requiring: 


*P  =  (C-10) 


The  n  numbers  are: 


q0t2 

ost2 

*2pc 

*3  - 

*  Pc 

*4  ~ 

*2Pc 

<M2 

at2 

fct2 

X2pc 

*6  ~  1 

*  Pc 

*7  ~ 

*2P, 
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K,t2 

Kjf 

*8  =  , 

*  Pc 

*9  = 

*2Pc 

u 

a 

*< 

II 

o 

Ik 

Ect2 

Ert2 

P* 

*n  =  2 

*  Pc 

*u  ~  , 

^  Pc 

U 

G. 

II 

H 

Ik 

vt 

at2 

gt2 

*15  = 

X 

*■16  =  — 

X 

*n  - 

X 

*0 

A, 

1  - 
II 

O 

Ik 

1  * 
n 

8 

Ik 

When  variables  with  the  same  dimension  are  scaled  in  the  same  way,  the  number 

of  scaled  variables  is  eight  and  the  number  of  the  non-dimensional  w  numbers  can 

be  reduced  to  five: 

at2  vt 

M  F 

pA2 

*1  —  »  *2  ~  » 

X  X 

*3  —  >  *4  —  » 

Aut2  X2a 

*5  - 

at2 

Table  C-l  summarizes  three  possible  prototype  -  model  ratios:  Each  column 
represents  a  different  set  of  possible  scaling  ratios  for  an  undistorted  system.  Note 
that  for  each  column  only  three  basic  ratios  can  be  determined  independently  (bold 
letters  in  Table  C-l).  The  ratios  listed  in  the  remaining  rows  are  determined 
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according  to  the  three  basic  independent  ratios  and  the  above  x  numbers.  According 
to  the  assumptions  listed  above,  the  geometrical  scaling  ratio  is  set  to  n,  and  the 
stresses  in  the  model  and  prototype  systems  are  the  same.  The  third  independent 
ratio  might  be,  for  example,  an  acceleration  ratio  of  one,  indicating  identical 
accelerations  for  the  model  and  the  prototype  (second  column  in  Table  C-l).  Thus 
proper  modelling  of  the  gravitational  acceleration  for  systems  that  do  not  involve  the 
means  to  change  it,  such  as  a  centrifuge,  is  ensured.  Note  that  in  this  case  the  mass 
densities  of  the  model  materials  should  be  increased  by  a  factor  of  n  in  order  to  keep 
the  scaled  model  undistorted.  If  the  mass  densities  are  kept  identical  (as  shown  in 
the  third  column  of  Table  C-l)  then  any  acceleration  of  the  model  system  should  be 
decreased  by  a  factor  of  n  when  extrapolated  to  a  prototype  system.  Similarly,  the 
prototype  time  is  larger  than  the  model  time  by  a  factor  of  n,  the  prototype  velocities 
are  equal  to  the  model  velocities,  and  any  prototype  force  should  be  extrapolated 
from  the  corresponding  model  value  by  increasing  it  by  a  factor  of  n2. 

Since  the  gravitational  acceleration  in  the  current  system  is  g= 386.4  in /sec2,  and 
at  the  same  time  the  mass  densities  of  the  soil  and  of  the  micro-concrete  are  similar 
to  the  prototype  materials  of  a  similar  larger  scale  system,  an  extrapolation  of  the 
experimental  results  to  a  larger  scale  system  involves  some  distortion.  However,  the 
dynamic  nature  of  this  system  makes  the  effect  of  gravitational  acceleration  relatively 
negligible,  and  reduces  the  distortion  in  the  range  of  other  experimental  errors. 
Evaluation  of  the  distortion  can  be  demonstrated  through  the  following  examples: 
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Example  1:  Effect  of  the  gravitational  acceleration  on  small  scale  model 
predictions  in  a  single  degree  of  freedom  system. 

Consider  a  single  degree  of  freedom  mass-spring  system.  Under  the  above 
assumptions  and  if  modelling  a  larger  scale  similar  system,  the  prototype-model  ratios 
are  Xp  =  n km,  and  gp  =  gm  =  g.  The  equation  of  motion  is: 

Ma  +  Kx  =  Mg  (C-ll) 

where  "a"  is  the  acceleration  of  the  mass,  V  is  the  displacement,  "M"  is  the 
lumped  mass,  and  "K"  is  the  spring  constant.  The  spring  constant  has  a  dimension 
of  Force/Length.  The  displacement  of  the  mass,  V,  can  be  written  as  a  summation 
of  two  components:  a  time  dependant  component,  and  a  static  component,  x*, 
which  depends  on  K,  M,  and  g: 

x  =  xdy  +  xst;  where:  x<,y  =  x^t),  x,,  =  x,t(K,M;g)  (C-12) 

Substituting  from  equation  C-ll,  yields: 

Mady  +  Kxdy  +  Kx*,  =  Mg;  x,t  =  Mg/K  (C-13) 


where  x,,y  =  f(o,t, initial  conditions,  and  o2  =  K/M. 
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It  is  assumed  that  the  similitude  requirements  following  the  dimensional  analysis 
for  the  soil-structure  system  also  hold  for  this  simple  system.  Hence,  for  a  true, 
undistorted  model  the  mass  densities  ratio  should  be  pp  =  pm/n,  which  follows  Mp 
=  n2Mffl  and  is  a  consequence  of  the  first  column  in  Table  C-l  that  can  be  checked 
as  follows:  Certifying  that  the  spring  constant  is  modelled  properly  (K  =  F/x, 
therefore  Kp  =  n2Fm/(nxm)  =  nK^),  the  prototype  displacement,  Xp,  can  be  written 
in  the  model  terms  as  follows: 

Xp  =  Xpdy  +  Xp*  =  nxmdy  +  Mpg/Kp  =  +  n2Mmg/(nKm)  =  n(xmdy  +  x^,) 

=  nxm  (C-14) 

The  distortion  is  caused  if  pp  =  pm  (or  Mp  =  n3Mm),  ap  =  a,„/n  but  "g"  remains 
unsealed.  From  x,,,  measured  in  the  experiment  it  follows: 

Xp  =  nxm  (C-15) 

True  scaling,  however,  yields: 

Xp  =  nxmdy  +  Xp,,  =  nxmdy  +  Mpg/Kp  =  nxmdy  +  n3Mmg/(nKn))  = 

nXmdy  +  n2Xnu.  =  +  fliflllKs.  (C-16) 

where  the  underlined  term  is  the  distortion,  which  results  from  the  model  to 


prototype  extrapolation  under  the  conditions  listed  above.  The  percentage  of  the 
distortion  can  be  defined  as  follows: 


%  Distortion 


•100  = 


1  + 


100 

1 

(n-1) 


(C-17) 


From  equation  C-17,  it  can  be  seen  that  if  the  static  portion  of  the  total 
displacement  in  the  model  is  less  than  1/200,  and  the  geometrical  scaling  factor  is 
less  than  50,  the  distortion  in  the  displacement  is  less  than  20%. 

Example  2:  Effect  of  the  gravitational  acceleration  on  small  scale  model 
predictions  of  the  displacement  of  the  roof  center  in  the  current  system. 

The  distortion  in  the  model  to  prototype  prediction  of  the  roof  center 
displacement  can  be  evaluated  by  assuming  a  dominant  first  mode  motion  of  the  roof 
slab.  In  this  case  equation  C-17  can  be  applied.  Consider  a  micro-concrete  circular 
roof  slab  embedded  in  the  sand  at  a  depth,  dm,  with  a  thickness  1^,  modulus  of 
elasticity  Em,  and  a  Poisson’s  ratio  vm.  Typical  values  of  these  properties  are  listed 
below,  and  yield  a  slab  stiffness,  Dm: 
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=  0.38" 

Em  -  2-106  psi 
vm  =  0.2 


-  Dm  =  Eh3/12(l-v2)  =  9830  lb-in 


(C-18) 


The  clear  span  of  the  slab  is  5  inches.  Assuming  a  uniformly  distributed  static 
pressure  within  the  slab  elastic  range,  the  static  deflection  at  the  center  of  the  slab 
is: 


*mst 


Pms.Rm4 

64Dm 


(5  +  v) 
(1  +  v) 

1 


;  simply  supported 
;  clamped 


(C-19) 


The  maximum  experimentally  measured  displacements  for  that  slab  were: 


0.04  inch  ;  simply  supported 

0.01  inch  ;  clamped 


(C-20) 


The  maximum  static  uniform  pressure  (ignoring  static  arching)  which  is  loading 
the  roof,  p^,  in  equation  C-19,  is  equal  to: 


P»,  -  Ymtl'm  +  =  (150/12^)0.38  +  (104/123)dm  (C-21) 


where  ymc  and  Ym*  3X6  the  weight  densities  of  the  micro-concrete  and  of  the  soil, 
respectively. 
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Substituting  p,,^  into  equation  C-19,  and  x^,  from  equation  C-19  and  x^10**1  from 
equation  C-20,  into  equation  C-17,  yields  the  distortion  for  the  above  conditions  and 
material  properties  as  a  function  of  the  depth  d,,,  and  the  geometric  scaling  factor, 
n,  as  a  parameter  and  is  shown  in  figure  C-2  and  C-3  for  a  clamped  and  simply 
supported  slab,  respectively.  As  can  be  seen  in  these  figures,  the  distortion  in  the 
prediction  of  the  roof  center  displacement  does  not  exceed  10%. 

Example  3:  Effect  of  the  gravitational  acceleration  on  small  scale  model 
predictions  of  the  vertical  pressure  over  the  embedded  roof  in  the  current  system. 

When  extrapolating  a  small  scale  model  stress,  pm,  to  a  prototype  stress  pp,  true 
modelling  requires  (Table  C-l): 

Pp  =  Pm  (C-22) 

However,  when  writing  the  pressure  as  a  summation  of  dynamic  and  static 
portions,  true  extrapolation  yields: 

Pp  ~  Ppdy  Ppst  —  Pmdy  +  Y ms(^m)  —  Pmdy  ^Pmst  — 

Pmdy  +  Pm*.  +  (n-l)Pm*.  =  Pm  +  (n-l)p„«  (C-23) 


Applying  the  previous  definition  for  the  percentage  of  distortion,  an  expression  for 
the  pressure  distortion  is  obtained,  similar  to  the  displacement  distortion  (equation 
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C-17): 


%  Distortion 


Pm+(n-\)Pma 


•100 


1  + 


100 

1 

(n-1) 


(C-24) 


Figures  C-4  and  C-5  describe  the  pressure  distortion  percentage  for  a  total  peak 
experimental  pressure  of  100  psi  and  25  psi,  respectively.  These  represent  some 
actual  peak  pressure  values  that  were  measured  in  the  tests.  The  static  portion  of 
the  pressure,  p^.,  was  expressed  as  a  function  of  the  depth  of  burial,  dm,  according 
to  equation  C-21.  It  can  be  seen  in  these  figures  that  in  the  range  of  depth  of  burial 
which  was  used  in  this  research  (2.5  to  3  inches)  the  distortion  in  the  peak  pressure 
predictions  is  less  than  30%,  for  a  geometrical  scaling  factor  as  high  as  60,  and  less 
than  20%  for  depth  of  burial  of  2.5  inches  and  a  scaling  factor  of  40. 

Using  a  similar  procedure,  it  can  be  shown  that  in  general,  the  model  to  prototype 
prediction  of  any  variable  in  the  system,  "VAR",  is  distorted  due  to  lack  of 
gravitational  acceleration  scaling  in  the  model  system,  according  to  the  static  portion 
of  that  variable  and  to  the  geometrical  scaling  factor,  n: 
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%  Distortion 


VAlC'-VAR 
— - MOO 

VAR^ 


100 

t  ,  1  VA  Rm 

+  (n-1)  VAR^ 


(C-25) 


Figures  C-6  (a  and  b)  illustrate  how  the  distortion  changes  with  the  geometrical 
scaling  factor,  n,  and  with  the  static  portion  of  the  variable  VAR,oul/VARstatlc. 
Figure  C-6(a)  describes  the  full  range  of  this  ratio,  and  figure  C-6(b)  the  typical 
range  in  the  current  research.  It  can  be  seen  that  the  distortion  due  to  gravitation 
decreases  as  the  dynamic  nature  of  the  phenomenon  increases,  and  as  the 
geometrical  scaling  factor  decreases.  It  can  also  be  seen,  that  the  distortion  in  the 
current  experimental  system  falls  well  within  other  experimental  errors,  and  therefore 
is,  indeed,  insignificant. 


m  =  i) 


229 


UOI}JO}SIQ  % 


o 
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ig.  C-2:  Distortion  in  displacement  modeling 
-0.38  inch,  x„,(total)=0.01  inch  (clamped  slab). 


Depth  of  Buri 

Fig.  C-3:  Distortion  in  di 
h=0.3B  inch,  xm(total)=0.04  inc 
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.  C-4:  Distortion  in  stress  modeling 
h=0.5  inch,  <Tm(total)=  100  psi. 


Depth  of  Burial  (inch) 

C-5:  Distortion  in  stress  modeling 
h=0.5  inch,  rrro(total)=  25  psi. 
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milliseconds 

Impulse  records  over  a  SIMPLY  SUPPORTED  roof  in  test  2  (4th  hit). 
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milliseconds 

Impulse  records  over  the  simply  supported  roof  in  test  6  (4thhit). 


D-7:  Pressure  over  the  clamped  roof  in  test  8  v.s.  free  field  (3rthit). 


Fig.  D-8:  Response  of  a  clamped  roof  in  test  8  (3rt  hit). 
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Impulse  records  over  a  clamped  roof  in  test  8  (3rd  hit). 
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